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1 Introduction 

The purpose of this work is to demonstrate that the noncommutative Kubo formula |]6] |53] |54] 21 |6T] |2l H] [12]| 
for the transport coefficients of aperiodic solids in magnetic fields can be efficiently evaluated on a computer. Our 
main message is that the noncommutative Kubo formula, which is now well known for its pivotal role in the 
theoretical foundation of charge transport in aperiodic systems ll53l [3] [4] |8] |2l |5] |40]|, can also provide an efi"ective 
numerical solution to the difficult problem of computing the transport coefficients for realistic quantum models of 
aperiodic systems. For disordered quantum lattice systems in magnetic fields, we derive a "canonical" finite-volume 
approximation to the noncommutative Kubo formula that can be efficiently evaluated on a computer. The advantage of 
this approximate formula is that it does not require twisted boundary conditions and integration over a Brillouin torus 
but only one (periodic-) boundary condition, and yet the convergence to its thermodynamic limit, which is the exact 
noncommutative Kubo formula, happens exponentially fast. 

The transport measurements have always been regarded as fine tools for probing and understanding the physical 
properties of solids but there have been particular times in solid state research when the transport measurements have 
taken the central stage. This was certainly the case when the Integer ||28l and Fractional [64| Quantum Hall Effects 
were discovered, or when the quasicrystals were discovered 155]. The recent discovery of the topological insulators 
ET\ l26l l27l [TOl [30l [39l [T8l U5 1 has brought again the transport measurements to the central stage. These materials 
are postulated to have non-vanishing transport coefficients even in the presence of strong disorder. A concentrated 
experimental effort towards proving this principle has generated an incredible amount of fine transport data for these 
materials. However, the goal has not been achieved yet, despite of several years of progress in sample preparation and 
characterization. As such, there is an urgent need to "reverse engineer" the available experimental transport data to get 
a better understanding of the microscopic structure of the samples. This is contingent on our ability to simulate the 
quantum transport for disordered materials in the presence of magnetic fields, based on microscopic quantum models. 
Even though we do not present any simulations for topological insulators, we wanted to bring the subject into the 
reader's attention because this was the primary motivation behind our effort. Simulations for topological insulators, 
using the formalism developed in this work, are reported in Ref. W^- 

The noncommutative Kubo formula for aperiodic systems appeared in Ref. fSJ, where the context was that of the 
Integer Quantum Hall Effect. It played a crucial role in explaining this effect. In the same reference, the authors derived 
the zero temperature limit of the Hall conductivity and linked it to the noncommutative Chern invariant. In the past 
several years, we have evaluated the noncommutative Chem invariant on a computer for several disordered quantum 
lattice models Il50ll59ll49ll48ll66l . These simulations revealed to us the extraordinary efficiency of this formula, which 
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Figure 1: The diagram summarizes the sequence of approximations. The random potential, represented here by the 
size-varying dots, is approximated first by a periodically repeating potential, and then the whole system is considered 
on a real-space torus. Comparative estimates are derived for the pair of algebras C'(Q x if, F) and C'(Qpgj. x Z'', F), 
and for the pair of algebras C*ia^^, x Z'', F) and C*iQ.A x T^, F). 



enabled us to compute this invariant for extremely large system sizes. We were also impressed with the accuracy of the 
approach for the quantization of the Chem invariant was observed to occur with better than five digits of precision (in 
the presence of strong disorder!), whenever the simulation box exceeded the localization length of the system. 

The results mentioned above stimulated our interest in the finite temperature noncommutative Kubo formula for 
charge transport. As it is always the case, on a computer we can only simulate finite systems. The noncommutative 
Kubo formula has earned many accolades, among which the fact that it explicitly tells us what the thermodynamic 
limit of the linear conductivity is. However, from a computational point of view, this formulation in the thermodynamic 
limit presents a major difficulty because many elements appearing in this formula do not have a canonical equivalent 
at finite volumes. This brings up a more general question, namely, how to accurately compute correlation functions for 
aperiodic systems on a finite volume? By correlation function we mean the expected values of products of observablesQ 

The solution to the above question can be constructed in a C*-algebraic setting. For this, we developed the 
concept of approximating C* -algebras, which are simplified versions of some "exact" C* -algebras, with the essential 
requirement that approximate morphisms can be established between the exact and the approximating structures. These 
approximate morphisms can be used to transfer the correlation functions from the exact algebra to the approximate one 
and to ultimately compute the error bounds for this approximation process. Our strategy will be to apply the concept of 
approximating algebras in two steps as illustrated in Fig.[T] The exact C -algebra is that of the observables defined in 
the presence of a uniform magnetic field (described by the antisymmetric tensor F) on the entire Z'^ lattice and on the 
entire disorder configuration space Q: C*(Q. x Z'', F). The first approximating algebra is that of the observables defined 
on the entire lattice but only on those disorder configurations leading to periodic potentials with a repeating cell A: 
C*(Qpef X Z'^, F). The differential calculus transfers automatically over C*(D.^^ x Z'', F) but a new trace is required. For 
the second step, C*(D.p^^ xZ'',F) becomes the exact algebra and the approximating algebra is that of the observables 
defined on a real-space torus of volume A and on a specific disorder configuration space Qf^: C*(Q\ x T^/, F). For the 
latter algebra we need to define an approximate differential calculus and a new trace. The approximate morphisms 
provide a canonical way to approximate the correlation functions, in particular, the non-commutative Kubo formula. 
For both approximating steps, we derive explicit error bounds on the finite-volume approximate Kubo formula and 
establish that the errors vanish exponentially fast as A is taken to infinity. This is the main result of our work. 

Since the results could be of interest to the condensed matter physics community and to the pure, applied and 
computational mathematics communities, we decided to make the presentation self-contained. The paper is organized 
as follows. In Section 2 we discuss the quantum lattice models for electrons in solids and we argue that, as long as 
one is interested in the linear response regime, these models should not be taken as some crude approximation of the 
physical reality, but rather as accurate representations of it, with a systematic way of improving the accuracy. In Section 

2 we also introduce an explicit model that will be used throughout the paper. 

In Sections 3 and 4 we review the noncommutative theory of charge transport in aperiodic solids developed by 
Bellissard, Schulz-Baldes and their collaborators |6, 53] |54l, touching only on the essential points. The interested 
reader can consult the comprehensive lecture notes by Bellissard 13 for a detailed exposition of the subject. Section 

3 introduces the noncommutative C* -algebra of the observables for a disordered electron gas in a magnetic field, 
and gives a short discussion of the functional calculus for C*-algebras. It also introduces the noncommutative 
differential and integral calculus on this C*-algebra. Together, these structures extend the notion of Brillouin torus 

*Since we are considering only the finite temperature regime, the observables in these products can be assumed to come from the analytic functional 
calculus with the Hamiltonian, which is assumed short-range. 
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to that of a noncommutative manifold, the famous noncommutative Brillouin torus flS). Section 4 reviews the basics 
of the noncommutative theory of electron transport and gives the derivation of the noncommutative Kubo formula. 
In both Section 3 and Section 4 we follow closely the work of Bellissard, Schulz-Baldes and their collaborators 



Section 5 develops the Combes-Thomas technique lfT4l for the noncommutative Brillouin torus. This technique, 
which is well known to the analysts working with the Schoedinger equation, will allow us to derive general and 
explicit exponential localization estimates for the analytic functional calculus with short-range observables. Section 6 
introduces the periodic C* -algebra C*(Qpgj. x Z'', F) and the noncommutative integro-differential calculus on it. It then 
presents the comparative estimates for the algebras C*(Q x Z'', F) and C*(Q.^^^ x Z'', F). 

Section 7 introduces the C*-algebra C*(Qa x Tj,F) over the torus. It also defines the noncommutative integral 
calculus and an approximate noncommutative differential calculus for this algebra. It then presents the comparative 
estimates for the algebras C*(D.p^^. x Z'', F) and C*{Q.a x T^, F). Section 8 introduces the approximate noncommutative 
Kubo formula on the torus and derives the error estimates, which are shown to vanish exponentially fast in the 
thermodynamic limit. Section 8 also discusses the numerical implementation of the approximate noncommutative 
Kubo formula. Section 9 presents a numerical application to the Integer Quantum Hall Effect and Section 10 
summarizes the main results of our work. 

2 Lattice models for solids 

The Hilbert space of a generic quantum lattice model is spanned by vectors of the form \n, a), where m e Z^ is a node of 
the lattice and a = 1, . . . , D is a label for the orbitals associated with that node. The scalar product (•, •) is by definition 
such that (n, a\m,/3} = S„mSa/3- Intuitively, one can think of \n, a) as the atomic orbitals of the elements present in the 
M-th repeating cell of a material. Formally, the vectors \n, a) can be viewed as the standard basis for the fi(Z'', C^) 
space. 

In the modern electronic structure theory, the quantum lattice models emerge as "low energy" approximations of 
the continuum quantum models. For perfectly periodic solids, there is a standard machinery for generating such lattice 
models from continuum models. In short, the construction starts from an effective Hamiltonian Hq which describes 
the equilibrium many-body state of the electrons in a material. A relevant spectral domain of Hq is identified, usually 
by defining an energy cut-off, and the subspace corresponding to this spectral domain is represented in the so called 
maximally localized Wannier basis functions (p„a 137 II . The number of Wannier functions per repeating cell depends 
on the number of atoms present in the repeating cell and on the energy cut-off. The vectors \n,a}, a = 1, . . . ,D, now 
represent the Wannier functions centered inside the n-th repeating cell and the lattice Hamiltonian is simply: 



where (•, •) is the scalar product of the continuum model. The above lattice Hamiltonian, although quadratic, it already 
incorporates interaction effects through Hq, which is tuned to exactly reproduce certain many-electron properties, such 
as the electron density for the case of the widely used Density Functional Theory [24 31]. The interaction effects, 
not accounted for already, are incorporated in the transport theory through the dissipation operator, to be discussed in 
Section 4. 

The individual terms inside the sums in Eq. [T]are called hoppings, and the matrix elements (0„,a,/^o<Am,/3) the 
hopping amplitudes. Generically, the maximally localized Wannier functions have a rapid spatial decay so the lattice 
Hamiltonians usually contain only a finite number of hopping terms between one lattice node and the rest of the 
nodes. It is quite often the case that just the hoppings between the nearest neighboring cells provide a reasonably 
accurate representation. If this is the case, we will use the terminology "nearest-neighbor hopping Hamiltonian." 
Another important observation is that the lattice Hamiltonian reproduces exactly the band spectrum of the continuum 
Hamiltonian Hq, below the energy cut-off. Furthermore, the coupling integrals with the external driving fields can 
be exactly computed for the Wannier representation and even a lattice representation of a realistic disorder can be 
generated |9 |. The main point of these observations is that we have a standard procedure to generate lattice models for 
periodic or disordered solids in external fields, whose accuracies can be systematically improved. Of course, the above 
statement is valid in the linear response regime where the representation based on the equilibrium Wannier functions 
makes sense. Several versions of this procedure are already implemented in the QUANTUM-ESPRESSO |20| first 
principle electronic structure calculations package, and applications to electron transport in molecular structures have 
already appeared in the literature 1331. 

The main conclusion of the present work, namely, that there is a canonical way to construct exponentially fast 
converging finite-volume approximations of the finite-temperature linear response coefficients, applies to a large class 
of disordered quantum lattice models. The essential characteristics of these models are: 1) existence of a triplet (a 
classical dynamical system) (Q, dPioj), t), where (Q., dP{oj)) is a compact metrizable probability space and t an ergodic. 
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probability preserving action of the group on Q, and 2) existence of a bounded family of lattice Hamiltonians /f„ 
indexed by the points of Q, with the property of being covariant: UaHi^U~^ = /^t„u- Here Ua represents the magnetic 
translation by a |6|. However, in order to build our final arguments, we will need to fully identify all the constants and 
the parameters appearing in our error bounds. For this reason, we decided to simplify the exposition by working with a 
specific lattice model. This model has just one orbital per node and, with disorder and magnetic field, the Hamiltonian 
of the model reads: 

neZ<'|pl = l neZ'' 

The phase factor e^^'" encodes the effect of the magnetic field through the Peierls substitution fA5\. Its specific form 
will be discussed in Section 3. Note that, in a short range isotropic model, all the hopping amplitudes must be equal in 
the absence of the magnetic field, in which case they can be set to unity like in Eq.|2]by choosing appropriate energy 
units. 

The last term in Eq. |2] 

V^ = wYj 0Jn\n){n\, (3) 

neZ'' 

represents the disorder potential. This is just a simplified but very popular representation of the real efi"ect of disorder 
in a solid. We will assume the amplitudes ol>„ to be independent random variables uniformly distributed within the 
interval [-5,5]- The collection of amplitudes {ca„}„eZ'' can be viewed as a point a> in the space O = [-1/2, 1/2]^'. It is 
now quite standard ||6l that Q is a compact metrizable space, which accepts the following probability measure: 

dPioj) = Y\ dco„. (4) 

neZ"* 

There is a natural action of the additive Z'' discrete group on D.: 

t„, : Q ^ D., (t,„oj) (5) 

which is ergodic and leaves the measure dP{co) invariant. Hence, the model of Eq.|2]gives a concrete representation of 
the generic quantum systems mentioned in the lines above. 



3 The noncommutative Brillouin torus 

In the Heisenberg's picture of Quantum Mechanics, the central objects are the observables, which are represented as 
elements of specific algebras. The main shift from the Schroedinger picture is the focus on the relations between the 
observables themselves, rather than on how the observables act on the wave functions. As such, the physics can be 
entirely formalized at a purely algebraic level. In this section, we follow Ref. ISl and introduce the C -algebra of the 
observables for a disordered electron gas in a uniform magnetic field, and the noncommutative differential and integral 
calculus over this algebra. We expedite the presentation and mention only what is strictly needed for the derivation of 
the noncommutative Kubo formula at finite temperatures. 

3.1 The C*- algebra for a disordered electron gas in a uniform magnetic field 

The Hamiltonian and its spatial translates generate an algebra which in some sense is the smallest (and therefore 
the most useful) space to carry the calculations in. It contains all the observables of interest and it provides just the 
right framework to compute the correlation functions at finite temperatures. Our first goal is to describe this algebra. 

The effect of a uniform magnetic field in a fi?-dimensional lattice model is captured entirely by a dxd 
antisymmetric matrix F [6|. For example, in 2-dimensions, \F\2\ is equal to the magnetic flux through the repeating cell. 
Our units will be chosen such that the magnetic flux is expressed in the units of magnetic flux quantum of (po = h/e. 
Also, we will use the notation jc • F • j as a shorthand for Yji,k^iP jkjk- 

The starting point of the construction is the algebra ^0, generated by the continuous functions with compact 
support defined on Q x Z'' with values in C, and with the following addition and multiplication rules: 

{af){oj, n) = af{a), n) (multiplication by scalars), 

(/ + g)(oj, n) = f{cD, n) + g{(D, n) (addition), (6) 
(/ *p g){cj, n) = 2^,2. fioj, m)git-Jaj, n - m)e"'^"" "'^ (multiplication). 

The index F attached to the multiplication operation is there to remind us that the multiplication depends explicitly on 
the magnetic field. The algebra JIq has an identity, denoted by 1 and defined by: 



l(w, n) = 6„fi. 



(7) 
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Every element / from Jio generates a covariant family of bounded linear operators {naif}a>en acting on £^('Z'): 

(Hz") 3cl,^ {{n^m in) = 2 f{t-'u, m - n)e"'^"'-P "'>cf>{m). (8) 

This is the representation of / e JIq as a linear operator acting on the quantum states, for a specific disorder 
configuration oj. The covariant property means: 



where Ua are the magnetic translations on ^^(Z°'): 
The map 7r„ is a (faithful) representation of J?lo algebra: 



(9) 



(10) 



(11) 



The representation of the identity on £^(Z'^), ;r^l, is simply the identity operator on £^(Z''). Any short-range 
Hamiltonian can be generated from an element residing inside this algebra. For example, the first-neighbor hopping 
model of Eq.|2] corresponds to the element: 



h(0J,n) = (5|„|,i(m) + WoinSnO, 

whose representation gives the covariant family of linear (self-adjoint) operators: 

neZ'i\p\=l neZ'' 



(12) 



(13) 



The Peierls phase factor can be identified now as (pp„ = n(p ■ F ■ n). 

The algebra Jio is too small to include interesting elements, such as the time evolution exp(-ith), which does not 
have a compact support. However, all the interesting elements can be generated from via limiting procedures and 
to rigorously define these limiting procedures we need to introduce first a norm. The norm on Jio is defined by: 



= sup ||7r„/||, 

(lien 



(14) 



where the norm on the right-hand side refers to the operator norm on {^(Z''y 



WQW = sup ||G«^||f2(z</), Il0llf2(z</) = 



2 \m\' 



1/2 



(15) 



The following is an upper bound on the norm of an element / e ^o, which is useful because involves only a simple 
estimate on the kernel of /: 



sup I V |/(w, m)| \ sup I V |/(t„w, n)\ 



1/2 



Quite often, we will use the following simplified version of this upper bound: 

11/11 < V sup|/(w,«)|. 



(16) 



(17) 



neZ'' 



aien 



Given that tt is a representation, and from the standard properties of the operator norm, the norm defined in Eq. 14 
can bee seen to have the following special property: 



ll/*F.?ll < 



(18) 



which makes (JIq, || ||) into a normed-algebra. Also note that ||1|| = 1. The completion of J?lo under this norm becomes 
a Banach algebra, and this Banach algebra is denoted by Ji. To make the Banach algebra Ji into a C*-algebra we need 
to define first a ^-operation, i.e. an operation with the following required properties: 



/"=/, (f + agy = f +ag\ (f*FgT=g**Ff, \\f\\ = 



(19) 
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Throughout our presentation, z will denote the complex conjugation of z e C. The following *-operation defined on J{ 



satisfies all the conditions stated in Eq. 19 



{n{co,n) = f{\-'oj,-n). (20) 
Furthermore, this * -operation satisfies the fundamental relation: 

ii/*frii = ii/ip, (21) 

which makes ^ into a C* -algebra. We collect all the above into the following central definition. 

Definition 3.1. The C* -algebra of the observables for the disordered electron gas in a uniform magnetic field 
(represented by the antisymmetric tensor F) is defined by the quadruple {J{,*f,\\ ||,*). We will use the notation 
C*(Q X Z'', F)for this C* -algebra. When there is no danger of confusion, we will just use 

The computations at finite temperatures can be carried entirely in the C* -algebra defined above. For zero 
temperature, one needs to extend algebra to a von Neumann algebra but since that is outside our scope, we stop 
here. 

Remarli. In the following we will drop the label F from *p and simply write the multiplication in J?l as / * 
3.2 Spectrum, the resolvent function, the analytic functional calculus and all that 

This subsection reviews some basic notions in operator algebras that are absolutely needed for our discussion. A good 
reference for this part is the short course on Spectral Theory by Arveson ||2J. 

Definition 3.2. The resolvent set of an element f belonging to an algebra with identity is the set: 

pif) = {z e C I / - zl is invertible}. (22) 

The spectrum off is the set: 

(^(f) = {z e C I / - zl « not invertible}. (23) 

Clearly p{f) = C - cr(/). 

The notion of the spectrum exists for any algebra with identity, but the characterization of the spectrum becomes 
more refined (and therefore more useful) if the algebra is endowed with additional structures. In particular: 

Proposition 3.3. For any element f in a Banach algebra with unit, cr{f) is a non-empty compact subset of the complex 
plane. Furthermore, cr(f) is located inside the circle of radius ||/|| and centered at the origin. If the algebra is a C*- 
algebra and f is a self-adjoint element, i.e. f* = f, then cr(f) is located on the real axis. Iff is a unitary element, i.e. 
f * f* — f* * f — 1' f'^^" 0"(/) located on the unit circle. 

The spectrum depends on the element itself but also on the algebra. As such, it is very common to use the wording 
"the spectrum of x in the algebra y," something that we will do quite often in the followings. Another central object for 
operator algebras is the resolvent function: 

Definition 3.4. Let f be an element of a Banach algebra J[. The following function: 

Pif) 3 z ^ (f - zir' e (24) 

is called the resolvent function off. 

Proposition 3.5. The resolvent function of any element f of a Banach algebra J[ is an algebra valued analytic function 
ofz e pif), in the sense that the limit 

,^(/-(z.01)-'-(/-zl)- (25) 

exists in J[ and is independent of how ( approaches the origin, for all z e pif)- 

Note that the above criterion of an analytic function, which we will use several times during our presentation, 
simply states that the function is differentiable and that the differential with respect to z is zero. Any other definition 
of an algebra valued analytic function is equivalent to the one considered above llT9l . For Banach algebras, the notion 
of complex integration over z of a continuously z-dependent element assumes the standard meaning via passing to 
the limit of the Riemann sum. Then complex integration and the resolvent function allows one to define the analytic 
functional calculus, i.e. a rigorous way to define functions of an element from the algebra. 
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Theorem 3.6. (The analytic functional calculus) Let f be an element of a Banach algebra !R and consider the algebra 
Tcr(f) of the analytic functions O in an open neighborhood of cr{f) (with the pointwise addition and multiplication). 
Then the following mapping: 

r^if) 3 ^ ^ ^ ^ - e ^ (26) 

is a homomorphism of algebras. Above, C is a contour surrounding cr(f) and contained in the domain of(b. The right 



hand side ofEq. 26 is regarded as <!>(/). 



The functional calculus can be pushed further to the algebra of continuous functions or even further to the algebra 
of bounded measurable functions. The latter will be needed if, for example, one does calculations at zero temperature 
where the Fermi-Dirac distribution has a jump at the Fermi energy. We however stop here, because we address only 
the finite temperature case. 



3.3 The noncommutative differential and integral calculus 

In the case of translationally invariant models, the transport coefficients are computed with the help of the integro- 
differential calculus with the ordinary functions defined over the Brillouin torus. When disorder and magnetic fields are 
present, the commutative algebra of the ordinary functions over the Brillouin torus is replaced by the noncommutative 
C*-algebra C*(0 xZ'',F). Remarkably, one can define a noncommutative integro-diff'erential calculus over C*(Q x 
Z'^, F), which transforms the C'-algebra into a noncommutative manifold, the famous noncommutative Brillouin torus 

m. 

For this, the A:-integration over the classical Brillouin torus is replaced by a trace over C*(Q. x Z"', F): 

T(f)= f dP(oj) f(a>,Q). (27) 
Jn 

By a trace, in general, is meant a linear functional with two basic properties: 

1. Positivity: T{f * /*) > for all / 0. 

2. Cyclic-city: 'r(/*g)='r(g*/). 



The definition of the trace in Eq. 27 is natural since, for the translational invariant case, T reduces to the ordinary 
A:-integration over the classical Brillouin torus. Another important observation is that, likewise the A:-integration over 
the classic Brillouin torus, the trace defined in Eq. 27 allows one to compute the correlation functions directly in the 
thermodynamic limit, in the sense that: 

lim ^JvMn^f)(n^g) . . .] = T{f * g . . (28) 

A->oo |A| 

where TrA means the trace over the quantum states inside the box A. The above relation is an immediate consequence 
of the Birkhofi" theorem and the ergodicity property of the translations on O |60|. Throughout our manuscript, \S\ will 
denote the cardinal of a set 5, and;^;^ will denote it's characteristic function. 

The derivations with respect to the ^/s of the ordinary functions defined over the Brillouin torus become 
automorphisms of algebra C*(Q x Z'', F): 

(djf)(u, n) = injfioj, n), j=l,...,d, (i= V^)- (29) 

Given a multi-index a = (ai, . . . ,ad), we will use the notation daf for d"^...d"j'f and we will write \a\ for 
a\ + . . . + ad. We collect some of the fundamental properties of the derivations in the following statement. 

Proposition 3.7. For j, k = I,. . .,d: 

1. The derivations commute: 



2. The derivations are *-derivations: 



3. The derivations satisfy the Leibniz rule: 



djdu = dkdj. (30) 
djif) = {djfr. (31) 



dj{f*g) = {djf)*g + f*{djg). 



(32) 
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4. The operator representation is given by: 

^Jdjf) = -i[xj,n^f], (33) 

where x = {x\, . . . , Xd) is the position operator on fiCif). 

It is important to note that the automorphisms 5/s cannot be extended to the whole C*-algebra C*(0 x Z'', B). As 
such, (9/s are not derivations in the strict C*-algebra language [52J. However, it is now customary to refer to (9/s as 
unbounded derivations (since they are un-bounded automorphisms in the true sense). Due to their simple action, the dj 
derivations can be integrated and one can easily find that they generate a c?-parameter group of isometric isomorphisms: 

\Xk-^^ (iii/)(w, n) = e-''' "f{co, n). (34) 

This group is not uniformly continuous or diflFerentiable of k over the entire The dense subspace of J[ containing the 
elements for which k uuf is A^-times continuously differentiable will be denoted by C'^{J[). An element / belongs 
to C^iJi) if and only if ||5q,/|| < oo, for all a with |Qr| < A^. Using the Leibniz rule and the fundamental property of the 
norm, one can easily show that if / and g belong to C'^(J?l), then their product f * g also belongs to C'^( Ji). In other 
words, C'^iJi) is a dense sub-algebra of €"{£1 x Z"', B). Finally, for any / e C\Jl), we have: 

dkMkf) = -djintf). (35) 

These technicalities will be needed when the quantum time evolution in the presence of an electric field will be 
discussed. 

A few important classical rules of calculus extend to the noncommutative calculus. 
Proposition 3.8. The following identities hold true: 
1. For f e C^(^) and f invertible in !K: 

djf-' =-f-'*{djf)*f-'. (36) 
2.If(S> and O' are any differentiable functions on the complex plane and f e C^{^), then: 

rmj) * dj<^'{f)) = 0. (37) 

This tells us that the noncommutative manifold has no boundaries, like the classic Brillouin torus f6^. 
3. (The noncommutative residue theorem l\47] . ) If<I> is an analytic function in a neighborhood of the unit circle and 
u is any unitary element, then: 

Tmu)djU) = a.iTiu-^djU), (38) 

where a-i is the coefficient ofz~^ in the Laurent expansion (S>(z) = Yj7=-x, <^nz"- 

It is also important for our discussion to mention that, once a trace is defined over a C* -algebra, one can define 
the LP norms: 

ii/iil/. = [T(|/|")]'^^ (39) 

and the corresponding L''{J{,T) Banach spaces, defined as the closure of J[ under the L'^-norms. A particularly 
important space is L^( Ji, T), since it can be endowed with a scalar product 

{f,8)=T(r*g), (40) 

which transforms L^(Ji,T) into a Hilbert space. As we shall see, certain automorphisms of the algebra Ji can be 
extended and viewed as self-adjoint operators on the Hilbert space L^(j?l,7"). An extremely fortunate consequence of 
all these is that we can use the functional calculus for self-adjoint operators on Hilbert spaces to define and compute 
functions of automorphisms. We will use this observation in our numerical analysis. 



4 The noncommutative theory of electron transport 

Here we give a telegraphic exposition of the electron transport theory in aperiodic systems, developed by Bellissard, 
Schulz-B aides and their collaborators ||6]|53]|54l. The goal is to present an explicit derivation of the noncommutative 



Kubo formula written in Eq. 87 in the most minimalist fashion possible. For a detailed exposition of the subject, the 
interested reader can consult the comprehensive lecture notes [5J by Bellissard. Before proceeding, we need a few 
more standard results in the theory of C* -algebras. 
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4. 1 Derivations and automorphisms 

Let g be a generic self-adjoint element of C*(Q x if , F), g* = g. Then g defines a bounded *-derivation (i.e. with the 
properties 2 and 3 of Proposition |3.7p via the following action: 

■Cgif) = i(8 * f - f * 8)- (41) 



The reverse statement is also true for a simple C* -algebra with unity 11521 . namely, any bounded *-derivation is 
generated by a self-adjoint element of the algebra via the action written above. Any bounded derivation generates 
a one-parameter uniformly continuous group of isometric automorphisms 1521 . To make this statement fully general, 
let us assume that g depends smoothly on the real parameter t and write g(t). 

Proposition 4.1. The equation for u( t, t'): 

dtu(t, f) = -X.g(t) ° u{t, t'), u{t' , t') = I ( = the identity automorphism) (42) 

has a unique solution, and this solution defines a one-parameter uniformly continuous group of isometric automor- 
phisms on C*{D. X 2f , F). In other words: 

1. u(t, t') o u(t', t") = u{t, t"), 

2. |Wf,0/ll = 11/11 

3. u(t, t) = 1, 

for any real t, t' and t". 

Below is another generic result in operator algebras, which will be important for our discussion. 

Proposition 4.2. Let u{t, t') be a one parameter group, invariant to the simultaneous translation of its arguments, i.e. 
u(t, t') = u(t - t'). If\\u(t)\\ < e"' for some a eM, then the Laplace transform of the one-parameter group is well defined 
and given by: 

f dt e-^t) = + £r\ (43) 
for all ^ e C with Re(^) > a. Above, X is the generator of the one-parameter group u. 



Proof. 

,->oo 

dt e-^'u(t) = lim V e-"^^'u(nAt)At 







n=0 



= |im^ ^ (e ^^'M(Af))" Af (44) 



n=0 

= lim Af(l - e-^^'u{At))-^ 

At^O 

= {i + £)-'. □ 

4.2 The quantum time evolution with and without a driving electric field 

As we already mentioned, the effective lattice Hamiltonians h considered in condensed matter physics involve a finite 
number of hopping terms from a site to its neighboring sites, thus is a self-adjoint element of J{q. We need to clarify 
that h describes the intrinsic properties of a material, i.e. it does not include the generators of the dynamics due to 
external fields. A short-range Hamiltonian defines a bounded * -derivation: 

Zh : C*(Q X Z^, F) C*(il x Z'', F), 

(45) 

Zh{f) = i{h*f-f*h). 

The one parameter group of automorphisms generated by X/, is the quantum time evolution in the absence of any 
external fields. It will be denoted by u{t, f) (or u(t - t')). Explicitly, u(t, t') is given by: 

C*(n xZ'^,F)3 f ^ u(t, t')f = * / * (46) 

where the elements e C*(Q. x Z"', F) are defined through the functional calculus. 

The quantum evolution in the presence of a uniform electric field E is not as straightforward. It is generated by 
the Hamiltonian: 

hE = h + Ex, (47) 
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and it will be denoted by ueQ, t') (or UE{t - t')). Here we set the electron charge to -1 (atomic units). This He is not 
an element of the algebra C*{Q. x Z"', F). Nevertheless, note that He still defines a well-behaved *-derivation over the 
algebra J{q: 

Zh,=Zh-EV,V = {du-..,dd), (48) 

which can be extended to the dense sub-algebra C^{J{) of C*{Q. x Z"', F). Furthermore, generates a one parameter 
group of isometric automorphisms on C*{Q. x Z^, F). This can be established as follows. 



Proposition 4.3. Let h e J?lo and h(t) = UieK with x\ defined in Eq. 34 Then h(t) is an element of for all t eR. As 

a consequence, the equation in v{t, f), 



dMt,t') = -£f,^,,ov(t,t'), v(t',t') = l. 



(49) 



(50) 



defines a uniformly continuous group of isometric automorphisms v(t, t') on C*(Q X Z'', F) (cf Proposition^^). Then: 

UE(t, t') : C*(Q X Z"', F) C*(Q x T^, F) 

UE(t, t') = U,£ o v(t, t') o UfE, 

defines a one parameter group of isometric automorphisms on C*{Q. X Z'^, F), whose generator is precisely X/i^. 
Proof. We start by showing that: 



■Ch(t) ='^tE° Zh°U,E- 



(51) 



Indeed: 



[v-tE ° £h ° U,e) if) = llr£ {■^■hiujf)) 

= iu,E [h * (u^j/) - (u-^f) * h) 

= iiutEh) * (u,£(it7£/)) - i{n,E{\x~^f)) * (u,£/j) 
= «((u,£/j) * f - f * {n,Eh)), 



(52) 



and the last line is precisely X.^t)f- Given the group properties of v(f, t'), it is straightforward to verify that ue defined 
in Eq.lSOlhas all the three properties listed in Proposition [411] Furthermore, we have: 



so the following equation holds true on C (Jl): 



UEit, t')-l= o v(f, f) o (UfE - U,e) + U,^ o (v(f, f') - 1) o U,e, 



UE{t,t')-l _i _i 

lim — j—p — = it,£ o (EV) o n,E - o 



(53) 



(54) 



If we use the identity from Eq. 5 1 and the fact that u commutes with V, we can see that the right hand side is exactly 
-(X/, - EV). a 



4.3 The collision processes 

Dissipation is a result of a perpetual sequence of electron scattering events. A scattering event can be described 
by a scattering potential w e C*(Q x Z'',F), whose exact form depends on the dissipation mechanism that is being 
considered. The scattering potential will, in general, fluctuate from one scattering event to another but for simplicity 
we will neglect such fluctuations. As discussed and exemplified in the lecture notes of Ref. (5\, w can be derived 
from microscopic many-electron quantum models describing the electron system, the environment (such as a phonon 
bath) and the coupling between the two system, by integrating out the degrees of freedom of the environment. General 
arguments show that w must commute with the Hamiltonian h. 

Following Refs. ||6l|53l|54l|5l, we introduce the collision processes via a time-dependent potential: 

w(t) = Y,S(t-tj)w. (55) 

A sequence of collision events can be indexed by a countable subset £ of R, where each point of £ represents the 
collision time of a particular scattering event from the sequence. The lengths sj of the time intervals between two 
consecutive collision events tj and tj+i are assumed to be independent random variables taking values in R+ = (0, oo) 
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and identically distributed according to the Poisson law e ^^''dsjT. t is called the collision time and it gives the average 
time between two consecutive collision events. Thus, the collision processes can be parametrized by a sequence: 

s = {...,s.i,so,si...). (56) 

The label will be attached to the time interval that contains the origin. To complete the parametrization, we need one 
more parameter, which is taken to be: 

g = dist{0,£nR+). (57) 

In other words, q is the time coordinate of the collision event occurring right after f = 0. As such, q takes values in the 
interval [0, sq]. With this parametrization, the collision times are given by tj = q + si + ... + sj-i. 
To summarize, the collision processes are parametrized by a point rj = (s,q) of the space 

H = (R^f X R+, (58) 

which is endowed with the probability measure: 



fifju(?7) = 



iX[o,saM)dql'd, (59) 



where X[o,sa]il) is the characteristic function of the interval [0, ^ol- We will use the shorthand dn{s) for e'^^'^ds/r 
and d/iis) for Yljex d^i(sj). A crucial observation is the existence of a natural action of the time translation at,t on £ 
which simply shifts each collision time by a Af. This action defines a one parameter group of automorphisms on the 
probability space S that are ergodic and measure preserving. As such, the time averaging and the ?7-averaging coincide 
||60 1 or, in other words, the collision processes are self-averaging. In the following, we will attach a label r] and write 



w,^{t) for the time dependent potential of Eq. 55 
4.4 The effective time evolution 

The time evolution of the system in the presence of a driving electric field and of the collisions events is generated by 
the Hamiltonian: 

h,^,E{t) = h + Ex + Wriit). (60) 

Since the collisions are assumed to take place instantaneously, we can write down the explicit expression of the time 
evolution: 



u,j,E(t, t') = UE{t- max{fj)) 



t>ti>t:-i>t' 



UEiminltj.i} -t'), (61) 



where the max/min is taken over the t/s appearing in the product. The automorphism w is defined as: 

wf = e"' * / * £>-"'". (62) 
Another way of writing the time evolution is (from here on we choose f =0): 

oo 

Ur,,E{t) = X(tifi)(t)UE{t) + ^X(tj^i,tj){t)UE{t - tj)w . . . UeUi - fl)H'M£(fi). (63) 

The time evolution is a covariant automorphism: 

U^,E(.t - Af, t' - At) = Ua^,^,E(t, f'), (64) 

and, as such, the averaging over the collision configurations: 

UE(t,t') = I driu„E(t,t'), (65) 



results in a one-parameter group of automorphisms that are invariant to time translations: UE(t, f) = UE(t - f). 

Our goal now is to compute the averaged time evolution M£(f). We will evaluate its Laplace transform, L[ue], 



starting from Eq. 63 Considering each term separately, we need to evaluate: 



I dte~^' I d^is) I dq/T X(tj,tj_i)(t)uE{t-tj)w...UE(t2-ti)wuEiti). 
Jo J Jo 



(66) 
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(67) 



This is same as: 

dnis) I dqlT I dt e~^'-'~''^UE(t - tj)w 
Jo Jtj 

X e-^^'^-'^-'^UE(tj - tj-i)w . . . e-^^'^-'''>UE{t2 - ti)we-^''uE(h) 
= J dfi(sj)J^ dte~^'uE(t)w J dfi(Sj-i)e~^'''-'uE(sj-i)w... 

J dfi{si)e~^'"uE{si)w J dfi{so) J dq/Te~^''uEiq). 

The emerging integrals can be evaluated explicitly, by directly applying Proposition |4.2| or slight variants of it: 

j dfx(sj) ' dt e-^'iiEit) = T (1 + + , (68) 

J dfi(sj) e-^'^UEisj) = (1 + T(^ + £„,))-', U = i-l,..., 1), (69) 



and 

^^0 



dq/T e-^^UEiq) = (1 + + X/,,))"'. (70) 

J Jo 
Thus: 

CO 

mw = T(i + T(^ + £h,)r' Yj ('^(i + + -Ch,))-')' , (71) 

7=0 



or 



L[u](0 = {^ + il-w)/T + £,y . (72) 



Given the statement of Proposition 4.2 we can now draw the important conclusion that the generator of the effective 
time evolution UE(t) is: 

(l-w)/T + £h,. (73) 
We introduce the notation F for the collision operator (1 -w)/t, and finally write: 

a£(f) = e-'(r+X*J^ ^74) 

The spectral properties of the collision operator have been discussed in Ref. JU for various dissipation mechanisms. 
In our numerical applications we will use a simplifying approximation for T, called the relaxation time approximation, 
where F becomes proportional to the identity automorphism. 

4.5 The noncommutative Kubo formula 

We assume that the electric field is turned on at f = 0, when the electron system was still in its finite temperature 
equilibrium state. In this state, the expected value (per unit volume) of an observable / e C*(0 x Z'',F) is given by 
T( f * po), where: 

po = (Dpd(/7) = (l + e'"-Ef^"^y' , (75) 

with <I>FD being the Fermi-Dirac distribution, T being the temperature, k the Boltzmann constant and Ee the Fermi 
energy. The Fermi energy is fixed by the electron density: = T(po). After the electric field was turned on, the 
expected values of the observables are given by T(f * p,), where p, = M;,,£(f, 0)po. As such, the charge current density 
at time t is: 

i,,£(0 =T(7*M;,,£(^0)po). (76) 



The charge current operator j is explicitly given by j = -Vh. Note that above we have a trace over volume, so Eq. 
gives indeed the current density. 

Now, the average over the collision events gives: 



76 



JeU) = J dfx{T]) T (j * M,,£(f, 0)po) = r (}• * UE{t)po) . (77) 
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The time average of this effective charge current density, 

Je = Km i f dt Jeit), (78) 
''^"^ t' Jo 

can be computed as: 

/-*oo 

Je = Hm 5 dt e-" T Q * S^COpo) - (79) 
''^o Jo 

At this step we see the Laplace transform of UE(t) appearing: 

/-*oo 

lim 6 I dt e'^' T (j * UE(t)po) =lim6T (j * L[u]{6)po) ■ (80) 



From Eq. 74 we can conclude at this step that: 

=lim5'r(7*(5 + r + X,J-'po)- (81) 
Since in the absence of an electric field the expected current is zero: 

5r(7*((5 + r + £„)"Vo) = o, (82) 

we can subtract such null contribution from Eq.[8T]to obtain: 

Je = * (,5 + r + X,,)"' o (EV) o(6 + r + Zi,T'pq) ■ (83) 

Finally, since Y and X/, commute with each other annihilate po, we see that the action of the last automorphism on po 
reduces to 6~^pa. At this point we can take the limit (5 — > to obtain: 

Je = -T ((Vh) * (F + £„y(EVpo)) . (84) 
The conductivity tensor o-jk is defined as the link between the charge current density and the electric field: 

d 

Jj = Yj '^J'^'" ^^^^ 
SO Eq. [84] already allows us to identify its expression: 

ajk(E) = -T ({dkh) * (F + Zh,y'dk^vD{h)) . (86) 

This expression depends on the applied field E through hence it is usually called the nonlinear conductivity tensor. 
The limit £ ^ defines the linear conductivity tensor which has the following expression: 

(Tjk = -T ({djh) * (F + X,)-'5,<I)fd(/j)) . (87) 
This is the famous noncommutative Kubo formula 161 15311541 . 

5 Analytic Functional Calculus: Exponential Localization 



We start here the discussion concerning the numerical evaluation of Eq. 87 The exponential localization estimates 
derived in this section are essential for the numerical implementation of the noncommutative Kubo formula. We will 
state our main results and the main conclusions for this section here at the beginning, filling the rest of the section with 
detailed proofs. 

Theorem 5.1. Let f e J?lo be an invertible element in !R. Then daf~^ is exponentially localized for any multi-index a. 
More precisely, let 

c(^)= sup V |(l-e-i«)|sup|/(^^,«)|, (88) 
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Definitely, c(^) is finite for any f e J?(o and c(^) is monotonically increasing of^. Let ^ > be defined by the equation 
c(f) = 1/11/"' ||. Then, as long as ^ < f , the following relation holds: 



\daf-\aj, n)\ < — ^, (89) 



where Cf^aiO <^ non-diverging parameter, in the sense that it takes a finite value when ^ = ^. In particular 

-m 

||/-i||-~c© 

and. 



ir'(^,«)i< „._i,i_i (90) 



|5,/-'(w, n)\ < '""^^^'"^ ^ (91) 

(ll/-'l|-'-c(^))' 

Remark. The parameter sup|j^(^j|^^ lk~'"^^'*^<9y/|| appearing in the last equation can be easily evaluated once / is 
explicitly given. For higher derivatives, explicit expressions for the coefficients Cf^aiO ^Iso available but they are 
relatively complicated. 



Corollary 5.2. Let h be the nearest-neighbor hopping Hamiltonian ofEq. 12 Then, for any multi-index a, there exists 
a finite, non-diverging parameter Ch,ai^), independent of z, such that: 

\daih - zIT\cj, n)\ < '^hAS'>e 

(dist(z,o-(/j))-2i/sinh^)l''l+i 



whenever 



In particular: 



and 



. . /dist(z^(A))\ 



e 



-m 



\{h - z\)-\u, n)\ < — . . ■ .. ^ (94) 

dist(z, cr(h)) - 2d sinh ^ 

^ (dist(z,cr(/z))-2t;sinh^)2 



Corollary 5.3. Let h be the nearest-neighbor Hamiltonian ofEq. 12 and let <P(z) be an analytic function in a complex 
neighborhood of cr{h) defined by dist(z, cr{h)) < k. Then there exists a finite parameter O, depending entirely on O, 
such that: 

|5„0(A)(a;, n)\ < ^'f/; , (96) 

(/f-2fifsinh^)l"l+' 

whenever 

^<| = s.nh-'(^). (97) 

Corollary 5.4. Let h and <l>(z) be like in the previous corollary. Then the operator representation n^<^{h) and its 
derivatives have exponentially localized kernels: 

\{nMa^{h))\m)\ < — , (98) 

for all £,<l. 

5 . 1 Thomas-Combes technique for the noncommutative Brillouin torus 



The main goal here is to prove Theorem 5.10 which is the key for exponential localization estimates stated above. Let 
^ be a c?-dimensional vector with complex components. We define the following map: 

e'f-^ : J?lo ^ (e'^^'fXoJ, n) = e'^"f(Lo, n). (99) 

Proposition 5.5. The map e'^'^ is an automorphism on which extends to the entire J{ algebra when Im(^) = 0. 
This automorphism sends 1 into 1. 
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(100) 



Proof. We have successively: 

m 

m 
m 

This confirms that: 

* g) = * (e^'^g)- (101) 

Clearly 

(e'^"l)(co, n) = 6„,o = Kto, n). (102) 

The map is invertible on Jio since e'^-^ o e"'^-^ = 1. Furthermore, if Im(^) = 0, then = ||/|| for all / e Jio, so 

e'f * can be extended by continuity over the entire Jl, as a norm-preserving map. □ 

Proposition 5.6. Fix f e JIq and consider e'^*/ as a function from with values in a normed space ( the algebra 
Then e'^'^f is analytic of^ on the entire C^. 

Proof. Consider a complex variation A^j of the j-th component of ^. This variation induces a variation of the 
vector ^ that will be denoted by A^j. We need to show that the limit 

hm ^ (103) 

exists and does not depend on how A^^- approaches the origin in the complex plane, for all ^ e C"'. This criterion for 
"analyticity" is equivalent to all the other existing criteria. We first guess that the Umit is equal to ixje'^'^f and then we 
prove that this is indeed the case. Note that ixje'^'^'f exists as an element of Jlo. We have successively: 

[(A^^)-i(e'(f^^^>V - e^-/) - ixje'^-'f] (co, n) 

= [(A^jr\e'"i^' - 1) - inj] e^''f{oj, n) (104) 
= 0(A^>^»/(w,«). 

Since / has a compact support, 0{A^j) can be bounded by an 0(A^^) that is independent of n. Then 

||(Af,)-i(e'<^^^^'-V - e^'V) - iXjc'^-'^fW < (5(A^,)|k^-»/ll ^ 0, (105) 
whenever \A^j\ is sent to zero. □ 

Proposition 5.7. If f e J{q is an invertible element in J{, then, for any real so is e'^^f. Furthermore, (e'^'*/)~' = 
e^-^f-\ 

Proof Even though f e f~^ will not be in However, since Im(^) = 0, the map e^'^ extendeds to the entire 
algebra J[, where it continues to be an automorphism. We are then allowed to consider the element e'^ */"^ Note that 
with the knowledge we have so far about this map, we will not be able to make sense of e'^'^f'^ if Im(^) ^ 0. Finally, 
we have: 

(e^*/)*(e'^V) = e^''C/*/"') = e^''l = l- □ (106) 

Proposition 5.8. Iff 6 J?lo is an invertible element in then there exists ^ > (to be given during the proof) such 
that e'^'^f is also invertible as long as |Im(^| < ^. 

Proof. Since 

and given the previous result, e'^'^f is invertible if is invertible. We have: 

^-im(f..)^ = (1 - (/ - * r') * L (108) 
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and, using Eq. [17 



IK/ - .'-■-(f-^V) * /-'II < lir'll sup V 1(1 - .--'-(f "')! sup \f{u, n)\. (109) 

Since the support of f{co, •) is contained in a bounded subset of "Z^ for all w's, the last line can be clearly seen to go to 
zero when |Im(^)| — > 0. As such, there is a f > such that 

ll(/-^'"'"<^"'/)*/"'ll<l, (110) 
for all ^'s with |Im(^)| < ^, in which case the sequence 

oo 

(111) 

converges in the alegebra J[ to an element that is precisely 

(l-(/-6'-i-(f->/) */-')"', (112) 
enabling us to conclude that e'^"f is invertible in these conditions. Furthermore: 

where 

ciO = sup V 1(1 - ^--i-^^"))! sup \f{a>, n)\, (114) 



is the constant appearing in Eq. 109 The upper limit ^ can be identified as the unique solution to the equation: 

c(i) = ii/-'r'. □ 

Proposition 5.9. Let us fix an invertible (in M) element f e Mq. Then the map: 

^-(^'f"/)"' (115) 

is analytic for all ^ e with |Im^| < f (with ^ defined in the previous Proposition). 

Proof. Considering again a variation A^j of the y-th component of ^, we need to show that the limit 

lim ^- ^ ^' (116) 

exists and does not depend on how A^j approaches the origin, for all ^ with |Imf | < f . We first guess that this limit is 
-(e'^'^ fyHixje'^'^ f)(e'^'^ f)~^ and then we prove that this is indeed the case. The inverse (e'^'/)"' exists for all ^ with 
|lm^| < |. We have: 

||(A^,)-'[(6''(f^^f'>-/)"' - (e'^'f)-'] + {e'^'f)-\ix,e'^-'fW^''fr'\\ 

= ||(^'(f-^f'>-/)-'6''f-[(A^,)-'(/ - e'^'^^^'^f) + ixjf](e'f-^fr' (117) 
- [(^''(f^^^''-^/)-' -(^''^■V)-'](/x,6''f-/)(6''f V)"'ll = 0(|A^,|), 



as we've seen in our previous estimates in Proposition |5]6] Thus, the limit in Eq. 116 exists and is equal to 



-{e'^''f)-\iXje'^\f){e'^\f)-\ □ (118) 

Theorem 5.10 (Thomas-Combes fT45 for algebra J{). Let f e J{q be an invertible element in Jl Then there exists a 
constant |, determined entirely by f, such that, for ally's with |Im(^) | < f , e^^'^f is invertible and we have the equality: 

=(e'f-/)-i. (119) 
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Proof. We have already learned that e'^'^ is an analytic automorphism on j?lo. Unfortunately is generally not 
in So we will force to be in J{q by applying the map;^^ : J{ — > J{q defined as (xRf){<jJ, n) = f{oj, n) if \n\ < R 
and zero otherwise. Then: 

^^^''^'CYfi/"'), (120) 

is an analytic map on the entire complex plane. Note that;^'ff commutes with e'^*: e'^'' o = Xr° ^'^^ on ■^o, and this 
property can be extended to the whole Ji by continuity. Note that we can always make sense of e'^^/"' as a complex 
valued function defined on H x Z''. As such, we can also say that: 

^^XRie'^'r'), (121) 



is an analytic map on the entire complex plane. Now, we have already seen in Propostion 5.9 that: 

^ ^ ie'^'fr' (122) 
is analytic for all ^ e C with |Im(^)| < ^. Obviously the map 

^^XRie'^'fr' (123) 
shares the same property. So at this point, we have two maps, 

XRie'^'f-') and XR(e'^'fr\ (124) 



which are analytic on the cylinder |Im(^)| < ^ and, according to Proposition 5.7 they coincide on the real axis. 
Therefore the two maps are identical on the entire cylinder |Im^| < f . Note that this statement is i ndepe ndent of how 
large R is. We can take the hmit 7? — > oo for XR(e'^^f)~^ because we already know (see Proposition 5.8 1 that {e'^^f)~^ 
exists as an element of Ji. In the light of the equality we just established, it follows thatXR(e'^^f~^) also has a limit in 
^ as 7? is taken to infinity, and this limit must be equal to (e'^^f)~^. □ 



5.2 Proof of Theorem ISTT] and its corollaries 



From the upper bounds on (e'^^f) ' established in Proposition 5.8 

-1 



we have 



or: 



(lir'ir' - c(^))" > Wie-^'fT'W = Ik^Vll > \e'^"f-'(oj,n)\, 

^Im(f «) 

\f-'(c,n)\< 



II/- 



■iii-i 



■c(^) 



(125) 



(126) 



independently of the direction of vector ^. Above we used the simple fact that |/(w, n)\ is always bounded by ||/||. Then 
the upper bound of Eq. 90 follows by taking Im(^) as -^j^, with ^ < f . 



To demonstrate the upper bound of Eq. 91 we make use of the identity from Eq. 36 

djf-' = -f-'*(djf)*f-'. 

Since e'^^ is a morphism, we can write: 



to conclude: 



W^"{djf-'){aj,n)\ < |k'f-^(5,/-')|| < |k'f-^5,/|| Ik'^V'll'- 
Using again the upper bound established in Proposition |5. 8 [ we can conclude: 

^ ^UP|j^(^)|^^ lie OjJW 



\(djf-')(uj,n)\< 



(ii/-'ii-^-c(f)r 



(127) 

(128) 
(129) 

(130) 



and Eq. 91 follows by taking Im(^) as -^j^, with ^ < ^■ 

Estimates on the higher order derivatives can be obtained in a similar way. For example, taking another derivative 



on Eq. 127 



- r' * (dkdjf) * /-' + * (djf) * * (5,/) * f-\ 



(131) 
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and following the same steps as before, we obtain: 

, 2||£'-i'"«-^>(5^/)|| lk-""<f-''(5^5,/)|| 

\e'^"{dkdjf-^){a), n)\ < ^ , + '■'!, . (132) 

{Wf-'W-'-ciOf {\\r\\-'-c(s)f 

This is the same as: 



|( W'Kw, n)\ < (133) 

(ll/-'l|-'-c(^))' 



with 



CkM) = sup (2||e-'-(f-)(5,/)|| \\e-"-^^-^\djf)\\ 

|Im(f)|=f 



+\\e 



(134) 



At this point we can take Im(^) as -^j^, with ^ < ^ and the desired estimate follows. The above parameter can be 
easily computed once the element / is explicitly given. Higher order derivatives can be treated in the same way. 

Corollary 5.2 is just a particular case of the main statement. For the case when f = h- zl, with h the nearest- 
neighbor (self-adjoint) Hamiltonian of Eq. 12 we can compute the parameters explicitly. Indeed c(^) can be computed 
from its very definition, to be: 

c(^) = 2c;sinh^. (135) 

It is important to remark that c(^) is entirely determine by h (so it is independent of z)- Also \\(h - zl)~' || = dist(2, o-(h)), 
a well known fact in the spectral theory. It is also straightforward to show that 

sup |k""<f '>5//i-zl)|| <2^^cosh^, (136) 

|Im(f)|=f 

which is again entirely determined by h. Corollary |5.2| then follows. 
For Corollary |5. 3 1 we have: 



(D(/j)=— {h-zV)-'<^{z)dz, (137) 

with Ck the contour defined by dist(z, cr(h)) = k. Using the estimates already established for the resolvent, we can 
automatically write: 

\y^'d„<^{h)\\<- — ^'^-"^^^ ^ r mz)\\dz\, (138) 

and the statement follows, with: 



O = ^ ^ mz)\\dz\. (139) 



Throughout the following sections, we will consistently use the notation from Eq. 139 



For Collorary 5.4 it is enough to observe that: 

\{nMaHh))\m)\ = da(t>(h)(t-'co, m - n). (140) 
Then the statement follow from the previous locaUzation estimates. 



6 Approximating algebras: First Round 

In this section we construct the algebra of observables for the more computationally manageable situations when the 
on-site potential is random only inside a large super-cell A and this super-cell is periodically repeated until the whole 
Z'^ is filled. We call this algebra the periodic algebra. Note that no condition is imposed on the magnetic field, so the 
whole system is not periodic. While the periodic algebra is not the end of our string of approximation, we want to 
point out that, especially in 2-dimensions, explicit analytic and computer assisted analysis is still possible at the level 
of this algebra. This is because one can chose a special (asymmetric-) gauge to represent the matrix F, so that the 
phases appearing in the multiplication law affect only one spatial direction. In this case the system becomes periodic 
in the remaining spatial direction and one can use the Bloch decomposition with respect to that direction. Hence, the 
2-dimensional Schrodinger equation can be reduced to a 1-dimensional Harper's like equation, which has been studied 
computationally by various techniques ifTSl [T6l l36l l63l ISTI |62 56 57]|32]|35][T71. As such, the rigorous error bounds 
reported in this section (see Theorem |6. 3 |l can still be of interest to numerical analysts. 
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6.1 The periodic C* -algebra 

Let tg . be the translation by one unit in the j-th direction and let A/' be a positive integer. We define the subset c Q 
of periodic oj's as: 

Q^, = e Q I t(2yv+iy,w = oj, j = I, . . . ,d}. (141) 

The set fl^^. is closed (hence compact) and translationally invariant, in the sense that t„Lo e fl^^, for all oj e Q^j. and 
n e H'- . An w from O^,. can be constructed by periodically repeating its components w„ for n in the unit cell defined 



per 

as: 



Kn--{-N,...,N\'. (142) 



We define the periodic C*-algebra CiO!^^^ x Z'',F) in the following way. The elements are continuous functions 
f(oj, n) from Q^j. x Z'' to C. Since Q^,. is invariant to translations, the multiplication law in Eq.|6] the *-operation in 



Eq. 20 and the operator representations in Eq. 8[ defined for algebra Jl, make sense when we restrict the /'s to Q^^. 
As such, we will keep *f , the operation in Eq. 2D]and the representation from Eq.|8]as the multiplication, *-operation 
and operator representation for the C*(Q^|. x '^\F), respectively. The norm is defined by: 

11/11 = sup ||7r„/||. (143) 

Obviously, ||/' * f\\ = \\ff, so C*(D.'^, x Z'',F) is a C*-algebra. 

The differential calculus can be defined as before because the automorphisms dj remain well defined when the 
functions are restricted to Q^j.. However, the the trace defined in Eq. 27 for algebra C*(Q x Z'',F) is of no use here 



because the set has zero measure relative to dP(a)). We define a new probability measure (only on fl^i), through 
the following natural formula: 

dPpAco) = Y] dcon. (144) 

Clearly this measure is invariant to the translations t because the w's are periodic, a property that allows us to define a 
trace over C*(Q^, x Z'^, F). 

Proposition 6.1. The following linear functional over C'(Q^f x Z'', F).- 

Tpe,.(/)= r fi?/'per(w)/(w,0) (145) 

has the following properties: 

Tper(/ * g) = Tperig * /),'rper(l) = 1, and TpM * /') > for / ^ 0. (146) 
As such, Tptr defines a trace over C*(Q^j. X Z'', F). 

Proof. The statement follows via a direct computation, using the fact that dPpeAc^) is translational invariant. □ 

Proposition 6.2. The following is an alternative, useful way to compute the trace: 



Tper(/)= f dPperico)Ti{in^f)XAj, (147) 

where ^^e characteristic function o/Ayv- 

Proof The statement follows from the definitions and from the fact that dPpet(a>) is translationally invariant. □ 

6.2 Analytic Functional Calculus: Comparative estimates 

Let fj e C*(D. x Z'', F), j = I, . . . ,N,he defined via the analytic calculus with the Hamiltonian: 

= 5„,(D//t), (148) 

with O/s analytic functions in a closed neighborhood of (r(h) defined by dist(z, cr(/j)) < k. Generally, we are interested 
in computing expected values of the form: 

f N 



T 



(149) 
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because most, if not all, the linear physical response coefficients of a system to external fields can be written in this 
way at finite temperatures. 

Now, any / from algebra J?l can be associated with an element from algebra C*(0^f x if, B) by just restricting its 
domain of definition from Q x Z'' to x Z"'. Let us denote this map by ''Pa^. Since the restriction eflFectively occurs 
in the Q space, the analytic calculus enjoys the following special property: 



(150) 



As already mentioned in the short introduction for the present section, the correlation functions in the periodic algebra 
may be evaluated on a computer. Then a natural question arises, namely, what errors should one expect if the elements 
from C*(D. x Z"', F) are mapped into elements of C*(Q^f x Z'', F) via ''Pa^, and the analytic calculus and the expected 
evaluation are performed inside the algebra C*(Q.p^^. x Z'', F) instead of C'(Q x Z'', F)7 This is the main question for 
the present section and the answer to this question is given below. 



Theorem 6.3. Let fj 's be defined as above and assume h is the nearest-neighbor Hamiltonian ofEq. 12 Then, for any 
< ^ < sinh" ^/(•/2fif), there exists a finite, completely identifiable (at least in the asymptotic limit of large N ' s) factor 
25i(f, {a}) that is independent off^j's, such that: 



T 



( 



/ N 



per 



W^aJj 

yj=i 



<23i(^,{a)) 



n 



<1>, 



N-'e-"^^^. 



(151) 



Proof. Since all 0//t)'s can be generates as 

o#) = ^ <i>j(z)(h-zir'dz, 

where Ck is defined by dist(z, cr(/j)) = k, we can focus first on /'s of the form: 

fj = d,^(h-zjir\ 



(152) 



(153) 



with Zi e Ck- In the following, we will use h- z instead of h - zl in order to simplify the notation. 

We partition the box A;^ in 3'' smaller boxes (assuming (2N + l)mod 3 = 0), in which case, at the center of A^/', 
there is a smaller box A;^ surrounded from all sides by the remaining 3^^ - 1 boxes of the same size. For any a> e Q., 
let PA;v the projector whose action is to take the restriction of to to Aaa and to repeat it periodically until a periodic 
to is defined over the entire Z''. Clearly pa„ o Va^ = Pa^v' Pa^v indeed a projector over Q. 

Our starting point is the following difference between the expected values computed on the inner box Aaa and for 
the operator representation: 



lA^I 



-Tr i 7T,„ 



( ^ 



-Tr^;r„ 



f N \ 



XAn 



(154) 



The first task will be to derive a bound on this difference. After simple manipulations, this diflFerence takes the following 
equivalent expression: 



Tr 



k=l 



V j<k 



i^OJ fk - ^VA^,Ojfk)^P 



\j>k 



Xa, 



(155) 



A basic identity for the calculus with the resolvent gives: 



dain^h - z) ' - da{n^,^^Ji - z) ' 
= da [{nji - zT\n^,i^^Ji - nMi^PAfjOjh - z)"') . 



(156) 



A key observation here is that the difference n^^^^h - n^h is equal to the difference between the onsite random 
potentials, so the previous equation reduces to: 



da ((TT^h - zr\V,,^^ - V^)in,,^^h - z)-') 



(157) 
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Note that the support of Vp^^M - V„ is contained in the complement = Tf - t^f^ oi the box A/^'. Furthermore, since 
we are considering only onsite disorder potentials, all derivatives of Vp^^w - V„ are identically zero so, by applying 
the Leibniz rule, the previous equation reduces to: 



(158) 



Thus the difference written in Eq. 154 is a finite sum of terms of the following generic form: 

A 1 



1 



-Tr { n„ 



(Vp, 



V J>k ) 



(159) 



where k = I,. . .,N, and: 

= d„.(h-z)-\ (160) 

with a~ + = Q"^. 

We exploit now the exponential localization of the analytic calculus. For some arbitrary product of /'s, we have: 



/ (/^-2t/sinh^)l"^l+'' 



(161) 



with Ch^OjiS) defined in Corollary 5.2 Hence we can automatically conclude that, for an arbitrary w e Q: 



V i<k ) 



A" 



(o), n) < 
(w, n) < 



{K-2d sinh^)!' 
C/i,<(^) 



n I 



(K-ld sinh^y 



C„,„/^) 



|ff,| + l 



(- 



-m 



(162) 



2t/sinh^)l«'l+' 



whenever ^ < sinh ^(K/2d). In the following, we will denote the product of the two coefficients appearing in front of 
e-fl"! by C(^): 

'^h.a-S^)Ch.al(^)nj^kCh,aj(^) 

<-■(?) = ., , ^« I I — ■ (lo3) 

(K-2d&mh^f^'^^r-^^"'^ 

Since Ci,^a-(.^)Ch,a+(.^) ^ Ch,at(.^), we can actually simplify the expression of C(^) to: 



C(^) = 



n^=iC,,,./^) 



(/.-2^sinh^)^^i-2r.,Kr 



(164) 



which is important because now C(^) is independent of k and is determined only by the initial data. 

Now, as we have already seen in Corollary |5.4[ the exponential localization of the elements implies the exponential 
localization of the matrix elements of their operator representation. Since {m\Vai - Vp^^oj\m} is non-zero only if m is 
located in the complement A^ of the large box Af^ and K/w|y^j - Vp,^^oj\m}\ < W, we can conclude: 



-Tr I n„ 



( N \ 

Ha 

V;=l / 



-Tr -^TT,, 



( N \ 
V7=l J 



I N 



£21-^1 aO-^ y J]e-'-^^"-"'^\{m\V^-V,,^Jm)\ 

J l^^l „eA„ meA' 



(165) 



f N 



< w 



W=l 



HeA« '«eA^ 



-2^\n-m\ 



where 2'^'' counts the number of generic terms, like those in Eq. 159 generated by the expansion of the original 
difference. Following the notation from Fig. [2] 



I ^ 2d oo 

i_ y y g-2(\n-m\ _ ^ V V V V g-2(y/(N-"j+l<)-+\y\- 



\An\ 



(166) 
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Figure 2: The summation over nz in Eq. 166 is done by partitioning is domains. There is one domain for each face 
of Ayv- The domain for the j-th face is shown as a shaded region. Then one fixes a distance k from m\.o Kj^ and let m 
walk parallel to the j-lh face, by varying the coordinate y. The step is then repeated for all allowed fe's. 



Due to the symmetry, each term in the sum over j will contribute with the same amount. Also, we will let y take values 
in the whole Z''"' and use the inequality '\l2(a^ + b^) >a + b,\.o continue: 



Id ■sr~\ v~i 
~ |A.,| 2-t 2-t 2-t 



-V2f(AA-«|+^'+|7|) 



At this point we combine the sums over k and over y into one sum over a 7 in Z"', to continue as: 



d 



< ^ ^ ^-V2f(A/-«i+iri) 



6d 



IN 

3d 



— y 
+ 1 ^ 



N/3 

e ■■" ) e 

ni=-N/3 yen' 
2V2, 



V2f|7| 



(167) 



(168) 



where the last estimate holds in the asymptotic limit of large A/"s. Here Aj(^) is defined by: 

AM) = 2 

Using a basic property of Arf(f), we can actually simplify the above upper bound, to conclude: 



^ ^ ^ g-2£\n-m\ < 



3d 



(V2^)yv- 



_2V2 



Collecting everything together: 



(169) 



(170) 



4 



-Tr-^TT,, 



1 



( N 



\An\ 



-TrU 



( N \ 



^A. 



1 -i^N 



(171) 



At this point we have an exponentially small upper bound on the difference written in Eq. 154 but there is one 
unsatisfying aspect, related to the fact that the periodically repeating cell is Af^ while the trace in Eq. 154 is taken only 
over the smaller box Af^. We will correct this by using a special shuffling map s which, in each repeating cell, switches 
the inner box with one of the remaining 3'' - 1 boxes inside that cell. Such shuffling map induces a natural map on Q 
that leaves invariant the probability measure dP(cj). We will use the same symbol s for this map. It is important to note 
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that s commutes with the projector pa^, which is true because the shuffling occurs only between the small boxes from 
the same repeating cell. The shuffling map also induces a canonical map on i^CZ^), to be denoted by s. Then we have 
successively: 



J dP{u)i:r[{n, 



t/PMTr (tt, 



f)X 



(172) 



The conclusion is that, when averaging over oj, the above trace does not change its value when A;^ is shuffled for 
another small box, so we can use 3'' - 1 shuffling operation to patch the entire Aaa and to conclude: 



dP{aj)Tr\i7T^^^f)X, 



(173) 



Therefore: 



3d, 



-W 



( N \ 

2 21- 



Arf^i(^)C(^)yV-'e-'^^ > 

N 



dPioS) 



1 



( N \ 



-Tr-^;?r„ 



rTr { TTn 



T 



lA^I 



( N \ 



( N \ 
V;=l J 



(174) 



fj 

V;=l ) 



Xf^N 



Finally, the integrand of the remaining integral over o) is independent of the values of co outside the box K^i, so those 
components of co can be integrated out and what is left is just an integral over [- ^j''^'^' with the measure written in 
Eq. 144 Then, from Proposition 6.2 it follows that the remaining integral is precisely Tperdljii ^^^a«/j)- 



It /}'s are of the form given liTEq. 148 we have to perform one additional operation, namely, we have to apply 
the following operation: 



(175) 



on the difference we just estimated in Eq. 174 As before, Ck is the contour defined by dist(z, o"(/i)) = k. The upper 
bound derived in Eq. 174 is independent of the actual position of the z/s on the contour Ck, so we can automatically 
conclude: 



I N 



T 



f N 



I N 



per 



4 



]~[5„^(I)/'!Pa^) 

f N 



(176) 



A^^i(^)C(^) 



and this concludes the proof. The expression of the factor 25i(f , [a]), in the asymptotic limit of large N, can be easily 



read from Eq. 176 



7 Approximating algebras: Second Round 

Here we go one step further and wrap the repeating cell K/^ oi the periodic algebra into a torus. We define a new 
C* -algebra over a cf-dimensional discrete torus and establish approximate morphisms between the periodic algebra and 
the torus algebra. The main result of the Section is that the analytic functional calculus with short-range elements from 
the periodic algebra is well approximated by the analytic functional calculus with elements from the algebra over the 
discrete torus. 
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7.1 A C* -algebra over the discrete torus 

We call a discrete circle S\j a discrete sequence of equally spaced points on a circle in a 2-dimensional plane. The d- 
dimensional discrete torus will be (S]j)^''. We will use letters like p, q, etc., to specify points on the discrete torus, 
while reserving n and m exclusively for points of Z''. On this torus we pick an arbitrary point and call it the origin o. 
Given a point p of the torus, we can use a succession of rigid rotations along the S\j circles, to rotate the torus until the 
origin reaches the position of the point p before the rotations. After such action, each original point of the torus has 
been replaced by a different point, and this substitution defines a map Xp on the torus. It is clear that Xp{o) = p, and that 
Xp is independent of the specific sequence (which is not unique) of S^^ rotations used to bring the origin at position p. 
Furthermore, the rotations satisfy the following commutative group relation: 

° = = XpOXg= x,^q. (177) 

It is convenient to require that each S^^ circles contain 2N + 1 points, though this requirement is not essential. We 
introduce a coordinate system on by naturally un-roUing the torus onto the points {-N, . . . , 0, . . . , N}'' of Z"', with 
the origin pinned at 0. The coordinates of a point p eT^ will be denoted by Hp. The set {-N, . . . , 0, . . . , N}"^ c Z'' will 
be denoted by A;^. This A/^' will later be connected with the repeating cell of the periodic algebra, hence the notation. 

We let to — [cL)p}p^j\i be a sequence of identical independent random variables with a uniform distribution in 
[-1/2,1/2]. Then co can be viewed as a point of the probability space = [-1/2,1/2]''^''' endowed with the 
probability measure: 

dPjioj) = Y] dojp. (178) 

The rotations r introduced in the previous paragraph induce a group of automorphisms on Q/^, which leave the measure 
dPi(cj) invariant. We will use the same notation r for the elements of this group. 

The C*-algebra C*{Q.jm xTii,F) over the torus is defined as follows. The elements are continous functions 
/ : Q.N X Td ^ C and the law of composition is: 

(/ * ~g){u,p) = 2 f{u, qrg{x-'u, x-'p)e>''^"'-^ "^\ (179) 

The operator representations on ^^(T^) can be defined as: 

((.^J)<p) (P) = Yj /^^p'^' (180) 

Proposition 7.1. The operator representation n defines an injective algebra morphism: 

nM*~8) = i^m~8) (181) 
if and only if all the entries ofF are quantized as: 

2 

Fij = X integer. (182) 

2N+1 ^ 

Remarlt. This is in line with Zak's finding that the magnetic translations accept finite representations only if the above 
quantization is satisfied ||69]| . 

Proof. A direct calculation will show that: 

((nM*m)(.P) = £ f\x-;co,x;^qrg(x-'uX,'q>'"'"^^^^^^^^^^ (183) 

9.9' eT,, 

and 

((nJ)(na,m)(.P)= Yj />;'w>^;'«)§(r^''^.^9'9')^''"*"'''"''""''-^"'V(9'). (184) 

9,9' eTj 

The equality between the two expressions will hold if and only if: 

n^^iq' ■ F ■ n^-iq + Hq' ■ F ■ lip = riq ■ F ■ Hp + riq' ■ F ■ riq + 2 X integer. (185) 

We have: 

ny-^>q = nq-rip + cr, n^-i q, = riq, - rip + a' , (186) 
where a and o"' are vectors with entries equal to or +{2N + 1). Eq. 185 then becomes: 

cr ■ F ■ {lip - iiq,) + cr' ■ F ■ {iiq - Hp) - cr ■ F ■ cr' = 2 X integer, (187) 



which is clearly satisfied, for all allowed iip, iiq and only if F,/s are quantized as in Eq. 182 The injective property 
follows directly from the definition of n. □ 
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From now on, we will assume the quantization of F stated in Eq.|182[ We define the norm on C*{^m x F) as: 



= sup 11^^/11, (188) 



and, given the result proven in Proposition 7.1 this norm has the fundamental property: 

ll/*§ll< 11/11 lllll, (189) 
transforming C*(Q^' x T^;, F) into a Banach algebra. Furthermore, a ^-operation can be defined as: 



r{a>,p)= f(x-'oj,x-'o\ (190) 
and the n representation and the *-operation defined above satisfy the essential relation: 

7^ J' = i^jT- (191) 



As such, the norm defined in Eq. 188 has the fundamental property: 

ll/*/*ll = ll/¥, (192) 

which makes C*{Qu x Tj, F) into a C*-algebra. 

The magnetic rotations on €^(Jd) can be defined as: 

([/,0)(/7) = 6'"'<"'-^«''V(r^V)- (193) 



With the quantization condition of Eq. 1 82 the magnetic rotations give a projective representation of the rotation group 



r. Furthermore, the following covariance property holds: 

Uq{nJ)U-' = TT,,^/. (194) 
Lastly, we introduce a trace on C*{Q.f^ x Tj, F) by: 

TT(.f)= f dPtioj) f{oj,o). (195) 



Proposition 7.2. The linear functional defined in Eq. 195 satisfies all the required properties of a trace. 



Proof. It follows immediately from the invariance of the measure dPj(cL>) relative to the rotations r. □ 
Proposition 7.3. The trace can be computed via the equivalent formula: 

Tt(/) = ^ r dP^iu) Tv{nJ}. (196) 

Proof. It follows immediately from the invariance of the measure dPficS) relative to the rotations r. □ 

7.2 A bridge between algebras 

We will define two natural maps which will be uses to navigate between the algebras C'(n^|. x 2f ,F) and C*(Q.n x 
Tj, F). First, let us note that Q^^. and Q.j^ can be identified in a canonical way. From now on, they will be considered 
one and same. The maps are: 

P : C*(Q^{, X T^, F) Ci^M x T^, F), 

(197) 

(Pf){a),p)= f{oj,np) 

and 

/ : ai^lN X Trf, F) ^ C'(n^, X If, F), 

(JfKcO, n)=Y, Sn,nJ(0J, q). ^^'^^^ 

qeT' 

It is important to note that (Jf)(a>, n) is null for n outside of A^v- These maps have several other important properties, 
which we state in the following two Propositions. 



26 



PRODAN 



Ti 




Figure 3: A point p with Hp e SAyy ^nd the set for Ti. In this simple case, Y.p consists of only one point, the point 
q. The actions of r"^ onp and Hp are illustrated. 



Proposition 7.4. Tlie maps P and J satisfy the relations: 



Po J - lc"(n«xTrf,F) 



(199) 



Proof. Indeed: 

HP o J)f)ico, p) = {Jf){ca, Hp) = 2] dn^^njioj, q) = Roj, p). (200) 

qeZ'' 



For the second identity: 



„ f f(to, n) if « e Ayv 

((/ o P)f){cj, n)=\ 5„,„^(Jf)(aj, q) = (201) 

0, otherwise. 



The following rule of calculus will prove to be extremely useful. 

Proposition 7.5. Let f e C*(Q.u x Tj, B) such that f(a>,p) = Oifnp> If i.e. if p is not a first neighbor of the origin). 
Then: 

if * Pg)(aj, p) = ((//) * g)ico. Hp) + % (202) 

where 

I if rip i dhjM 
2 fioj, q)[g{i-"^oj, - g{t-"^co. Hp - ng)Y-("^""^\ if Hp e 5A^. ^^03) 

Here 'Lp is the set of those q's who are first neighbors of the origin and their corresponding rotations r~'/> send tip, 
which is necessarily located on one or more faces of hj^i, to an opposing face of (s^^ Fig- 

Proof 

{f^Pg){u,p) = J] f{aj,q)iPg){r-'co,x-'p)e-''^"'-""^^ 
= Y,fi^,q)8(i-"^oj,n,-^,p)e"^'-"^""^^ 

qeTj 

= Yj((P° J)f)('^' q)8(i~"'<^, n,-.p)e"'^"^" "^^ 
= 2(-^/)(^' ng)g(t-"^(j, n,-^,p)e'<"'-^ "->^ 

If lip ? (9Ajv, then = rip - tiq because riq is a first neighbor of and we can continue as: 

. . . = Y}Jf){cj, nq)g{r"^co. Hp - «,)6'-<«''-^«') 

= Yj (JDi^' m)g{r"'to, lip - m)e"'^"''" '"^ (205) 

meZ'' 

= ({Jf)*g){aj,np). 



(204) 
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If Hp e dA^, then the rotation can send n^-ip to the opposite facet of A^^, relative to Hp. Thus, the product is no 
longer equal to ((//) * gXo), Hp), but instead: 

(/ * Pg)ic, p) = ((//) * g) (oj. Hp) + 2 (Jf)(aj, n,)[g(t-"^cj, - g(t-"^oj, lip - «,)]e"^(«''-^«'>. □ (206) 
Remark. The above result can be used to demonstrate that / and P are approximate homomorphisms. 



7.3 The analytic functional calculus: Comparative estimates 



Let h e C*(D.x 11', F) be the Hamiltonian of a quantum system. Our string of thoughts was to restrict h to periodic 
(jj's via the map ^a„, and then to approximate ""^^Ji on the torus by ^ = Pi^i^^h). If there are only a finite number 
of hoppings from one site to the neighboring sites, then ^Afjh can be recovered from the Hamiltonian on the torus 
via ^Kfjh = Jh, so no truncation or loss of information occurs for the Hamiltonian itself during the last step of 
approximation. However, the functional calculus with the Hamiltonian will be difi"erent for the two algebras. Indeed, 
although ^^Afjh and h are similar as complex functions, the inverse of ^ - zl, which generates the analytic functional 
calculus on C*(Q.f^ x Tj, F), will be sought inside C*(Q.f^ x Tj, F), while the inverse of '^A^jh - zl, which generates 



N 

per ■ 



l'^, F), will be sought inside C*(Q^, x Z'', F). 



the analytic functional calculus on C'(Q 

Our task now is to demonstrate that the two generators of the analytic functional calculi converge to each other 
exponentially fast when the size of the torus is increased to infinity. To keep the technical complications to the 



minimum, we will consider only nearest-neighbor hopping Hamiltonian of Eq. 12 However, the main conclusion 



of this subsection, contained in Theorem 7.8 can be easily extended to any model with short-range Hamiltonian. 



Proposition 7.6. The following relations between the spectra always holds: 

cr(h) C cr( Jh). 



(207) 



Proof. With the quantization of Eq. 182 Jh becomes a periodic Hamiltonian on Z'^. Then one can decompose 
Jh according to the irreducible representations of the discrete translations, which are indexed by a A: point of the d- 
dimensional torus: h = Jh{k)dk. It is straightforward to show that h is equal to Jh(0). The statement follow because 
aiJh) = LlkO-[jh{k)). □ 



Lemma 7.7. Let h be the Hamiltonian on the torus corresponding to the nearest-neighbor hopping Hamiltonian of 
Eq. 12 and let z & C be outside the spectrum of Jh. Then: 



- zlT' - PiJh - zlT'W < 



N 



dist(z, cr(h)) (dist(z, cr( Jh)) - 2d sinh^) ' 

whenever^ < smh~^(dist(z, cr(Jh))/2d). 

Proof Let r e C^iD.^ x T^, F) be defined by 

l + f = (h-zl)* P(Jh - zlT^ 

According to Proposition [73] 

if «p ^ dAjsi 
r{to,p) = \ ~J] '^MJ'^ - zl)"kt""'w, w,^ip) 

- (Jh - zl)-'(t""''w, lip - «,)]£'"'<"''-^"'' if Hp e dAf^. 



(208) 



(209) 



(210) 



The number of points in is less or equ al to d , the equality occuring when rip is located at a comer of Aj^j. Note that 
the second argument of { Jh - zl)~' in Eq. 210 is on the boundary or outside A/^ and the key here is that {Jh - zl)~^ is 
exponentially localized. As such: 



\r{(^,P)\- 



= if rip i dAf^ 
2de-^\"-\ 



dist(z, o-{Jh)) - 2d sinh^ 



if rip e dA/^. 



(211) 
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To compute the norm, we proceed as follows: 

||r|| < ^ &upl\r(oj,p)\} 



< ^ y e-fi'vi 

dist(z,(r(//i))-2^sinh^„^^^ 



2:^ 2d y ^-fV^-+i'«i' 

-2cfsinh^ ^ 



(212) 



dist(z, cr(Jh)) - 2d sinh ^ 



dist(z, cr(//i)) - 2d sinh^ 



2 ^-4(A^+I".|)^ 



to conclude: 



4fi;%_i(^/V2)£'- 



V2 



N 



\\r\\ < " ' . (213) 

dist(z,o-(//!))-2cfsinh^ 



Since cr(k) c cr( Jh), we can multiply Eq. 209 with (h - zl) \ to obtain: 

(h - zl)-' - P(Jh - zl)-' = -(/i - zl)-' * r (214) 
and the statement follows. □ 

Theorem 7.8. The analytic functional calculi with the Hamiltonians corresponding to the nearest-neighbor hopping 



Hamiltonian of Eq. 12 in the C* -algebras C*{Ql^^ X 7f , F) and C*{Q.u X T^/, F), become exponentially close as N is 
taken to infinity. Specifically: 

mh) - P^mW < ' , (215) 

K(K-2d sinh^) 

for any function 0(z) analytic in the neighborhood dis\.{z, cr(Jh)) < k ofcr(Jh) and for ^ < sinh- ' (A-/2fif). 
Proof. The statement follows from the representation: 



= -[' 



<D(/i) - PO(Jh) = — I <D(z)[(/j - zl)-' - P{Jh - zl)-']fi;z, (216) 



and from the upper bound of Lemma 7.7 



7.4 An approximate differential calculus on the torus 

Our goal here is to demonstrate that an approximate differential calculus can be defined over the torus algebra such 
that, for elements generated by the analytic functional calculus with a short-range Hamiltonian, the errors introduced 
by the approximate derivations vanish exponentially fast as the size of the discrete torus is taken to infinity. 

Theorem 7.9. Let x : .Sjj M.be a continuous function such that (n^, is the coordinate of the point p e S\^): 

Hp) - np\ = if \np\ <N/^ (217) 

and, 

|i(p)| < Kl for all p e 51,. (218) 
We define an approximate differential calculus on the torus by: 

0jf){p) = ix{pj)f{p), j=l,...,d. (219) 

Here, pj denotes the j-th component of p. Let h be the Hamiltonian on the torus corresponding to the nearest-neighbor 



hopping Hamiltonian ofEq. 12 and let <I)(z) be an analytic function in the neighborhood dist(z, cr(Jh)) < k of cr{Jh). 
Then for any < ^ < sinh-'(A'/2c?), there exists a finite and fully identifiable parameter 932(^, a) such that: 

\\da(^{h) - Pda^{Jh)\\ < »2(^, a)5)V'"^'^ e'^-'^, (220) 
for any multi-index a = {a\,a2,--., ad}. 

Remark. It is very likely that the above estimates can be improved and the power law factor can be entirely eliminated. 
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Proof. We will use the notations: 

^aip) ■■= ^{pir . ■ ■ KPdT", and n,(p) := . . . (221) 

We have: 

[d„m)-Pd„<:)iJh)]{aj,p) 

= i^"kaip)<I>mco, p) - i^"^nJp)(PO(Jh))(oj, p) (222) 
= jl'^li„(/7)[(D(/t) - P<^{JmoJ,p) + i^"\^a{p) - na{p)mJh){cD, rip). 

From Theorem l7.8l and the conditions assumed on function i: 

\^a{p)[^{h) - P(^{JK)\{u,p)\ < ^, . • (223) 

' ' K{K-2a sinh ^) 



Furthermore, from Corollary [53] 

Xa^ (rip)^ \n„(p)\ ^--fl""! 

\M) - n„ip)mJh)(aj, np)\ < ■ ■ (224) 

' ' K -2d sinh ^ 

These estimates and the following coarse upper bound on the operator norm: 

\\da(^(h) - Pda<l>(Jh)\\ < {IN + lY sup |[5„<D(ft) - Pda<^(Jh)\(cD,p)\ , (225) 

allow us to conclude that, in the asymptotic limit of large N's: 

\\da<I>(h) - Pda<^(Jh)\\ < 2' \4d^Ad-i(^/ V2)/^ + 1 , ■ ■ . ^ , (226) 

L i K-2d smh ^ 

and the statement follows. The parameter 252(f, «) can be identified from the last inequality. □ 

It remains to determine how the approximate diflFerential calculus works in the operator representation and to 
explicitly construct and test a practical function i. 

Proposition 7.10. Let O be the discrete unit circle in the complex plane defined by the solutions of the equation 
= 1, and let 

x(p) = YjbAr^ (227) 
be the discrete Fourier decomposition of the function x. Then: 

~^Jdjf) = i 2] b,A-''(7:^f)A'', (228) 

AeO 

where Xj is the operator (Xj(f>)(p) = np.(p{p) on ^^(T^). 



Proof. Using the explicit expression of the operator representation given in Eq. 180 
[injdjfmip) = Yj^djf){x-/u,x-/q)e"'^"^-P"^^m 

= ^(('^;'«)7)/('--;'^'^-;'9)^"^"''^"''V(?)- 



(229) 



A key observation is that n^-ig = itg - Hp + cr where the components of cr are either or ±{2N + 1). As such, for all 
A eOwe can write: 

/('"''),= rT""^. (230) 
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Then we can continue: 



AeO 9eT</ 



(231) 



AeO 



YbA[r'^(nJ)(A''mp)- n 



AeO 



To construct a practical solution, it is convenient to perform a rescaling of the coordinate system on S^, (hence on 



Hp Xp := UplilN + 1). 
Note that the new variable Xp always take values in the interval [-^, i]. With this rescaling, 

where A takes the values e^'^"', with k = -N, . . . ,N. The function i needs to be rescaled too: 

x^x{p) = xip)/i2N+l). 
The Fourier series of the function x takes the form: 

N 



(232) 



(233) 



(234) 



(235) 



k=-N 



and since x{p) is real and can be chosen odd relative to the reflection relative to the origin, the Fourier series reduces 
to: 



1 ^ 

' k=l 



(236) 



A key observation about the above rescaling is that the A/^'s do not depend on N anymore (only the upper and 
lower limits of k depend on N). Furthermore, since the elements obtained from the analytic functional calculus are 
exponentially localized near the origin, the rescaling of the coordinate system concentrates these elements more and 
more towards the origin, as M is increased. Thus, a good approximation for the differential calculus is obtained by 
requiring i to generate a good representation of the function ix near the origin. In the same time, we need to make sure 
that the function closes continuously on the discrete circle. These two requirements are fulfilled by a Fourier series 
i Yjk=i ^kO-l - \'^)^ where the coeflicients bk are chosen such that the difference: 



1 ^ 

X - ^Yj^,ai-i7) - oix''-'). 



(237) 



k=\ 



asymptotically for x near the origin. Here, Q is taken to be a fraction of N. This condition is fulfilled if i>t's are 
solutions to the following system of linear equations: 



Q . 

Yk^J-^bt = -5j,,, je[\,...,Q]. 



(238) 



The matrix on the left side of the above linear system of equations is a Vandermonde matrix and, as such, its explicit 
inverse is known [34| and the linear system of equations can be easily solved. 
To summarize, we propose: 

Q 

x{p) = K2N +l)Y h{^7 - ^a""')' (239) 

k=l 
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20 40 60 



Figure 4: Left: Representation of the function x(p) = i{2N + 1) b^iA'l^'' - ^^i^"''), when Q is varied from 2 to 10 in 
increments of 2. Right: Plot of the difference Hp - x(p) for the same values of Q. 



with b^s satisfying the linear system of equation from Eq. 238 With this choice: 

xip) - Up ~ oiUpIN)^, 



(240) 



in the asymptotic limit rip « N, which is relevant for us. Furthermore x(p) closes smoothly on the discrete circle. 
Fig. |4] illustrates the behavior of x(p) as Q is increased from 2 to 10, for a typical system size (N = 60) used in 
our computer simulations. As one can see, with a good approximation, x(p) satisfies all the requirements stated in 



Proposition 7.10 when Q is given a reasonable value. In our numerical simulations we take Q = 10. 



7.5 Comparative estimates for the correlation functions 

Lemma 7.11. Let h be the Hamiltonian in C*(f2^f X Z'', F) corresponding to the nearest-neighbor hopping 



Hamiltonian of Eq. 12 and (j = I, . . . ,N) be a set of analytic functions in the neighborhood of cr(h) defined 
by dist(z, cr{h)) < k. Let fj (j = l,...,N) be elements of the C* -algebra C*(0^, X T^, F) of the form: 



fj=da.^j(h), j=l,...,N. 

Then for any < ^ < sinh~'(/f/2fi?), there exists a finite and fully identifiable parameter 233(^, [a]) such that: 



(241) 



^ N \ 



per 



t N 



V>=1 J 



Proof. In general: 



so we can write: 



r N \ 



Tr 



per 



f N \ 



( N 



- Tt 



N \ 



(242) 



(243) 



(244) 



We see that the problems has been reduced to comparing P Yl%i fj ^nd n^=i Pfj- We proceed as follows. Let r^ denote 
the difference: 

k f k-\ ^ 

rk = PWfj- PWf^ *Pfk- (245) 
j=l V j=l ) 



Then: 



N N N ( 

j=l j=l k=l \j>k 



(246) 



We show in the following that each r^'s are exponentially small. From definitions: 

fk-i \ (k-\ \ 



(u,ng)fk{t„laj, n,-^,p)e 



m(npFn„) 



(247) 
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Since n^-ip = rip - Uq if Hp - Hq e Af^, many terms cancel in the above difference, and all that remains is: 

fk-i 



rk{co,p)= Yj Yl-^J 

meZ'',«p-meA^ \j=l J 
(k-l \ 

- Z \Uf. 



{co,m)Mt-Jco,np-m)e""'"'-"-"^ 



(bJ, ng)Mt„^^aj,n,-^ip)e 



We now make use of the following bounds: 

fk-i \ 



(cD, n) 



k-l 



< e 



(k-2d sinh^)l«^l+^ 



and 



to write: 



\rk((^,p)\ < „ , „ X 

()t-2t;sinh^)2.^i<l"^l+^) 



-^(\m\+\np-m\ 



Z 



-f(|n,|+|n, pi 



We use the fact that the exponent \m\ + [rtp - m\ is larger than N, \m\ or to write: 

\m\ + \np -m\>^N+ ^\m\ + 



Similarly, + \n^^^p\ is larger than A/", |«^| or so 



1 . . 1. . 1. 



\ng\ + |Mv„pl > + -^\'^p\- 



As such: 



( k 

n 



()t-2fifsinh^) 



We can now estimate the norm of using the bound of Eq. 17 



Also: 



()t-2t/sinh^)2t,(W+i) 
Hp/;!! < C,,„,(0)5),7^l«'l^i. 



Finally, based on Eq. 246 



we can write: 



' N N " 



<Yj\\rk\\Y\\\Pfj\l 

k=i j>k 



and conclude that: 



is bounded by: 



Tr 



per 



' N \ r N 

\j=i J \}=\ J 



\k=\ i<k ^" "".^iii^*^/ .^1^ 

The parameter 233 W) can be explicitly read from above. 



y n ^^^^ n^^lfno,^ 

f-^i \ (/^-2t;sinhOl«''l+i 1 1 /fW+i 1 1 ^ 

A-l kA: ^ ;>^' / V ;=1 
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Theorem 7.12. Let h be the Hamiltonian on the torus corresponding to the nearest-neighbor hopping Hamiltonian of 
Eq. 12 and <I>j (j = I, ... ,N), be a set of analytic functions in the neighborhood ofcr(Jh) defined by dist(z, cr(Jh)) < k. 
Let fj,j= I,. . . ,N, be elements of the C* -algebra Ciil^ X T^, F) of the fonn: 



fj = ~da.^j(h), j=l,...,N, 

and let fj,j= I, . . . ,N, be elements of the C -algebra 

C*(n^r X Z'', F) defined by: 
fj=d„^<i>jiJh), j=l,...,N. 
Then for any < ^ < sinh~'(/<'/2fi?), there exists a finite and fully identifiable parameter 234(^, [a]) such that: 



per 



{ N \ ( N \ 

Uf^ -^T 

yj=i J \j=i J 



< 254(^,{a)) 



f N \ 



n 



e 2 



N 



Proof. We write: 



per 



( N > 

\j=l J 


- Tt 


( ^ 




( ^ ^ 

1 7=1 ) 


- Tt 


( 

b=i J 


V 



( N ^ \ 

Vj=l ,7=1 ) 



(260) 



(261) 



(262) 



(263) 



The absolute value of the first difference, |Tper(n^=i fj) " '^T(n^=i Pfj)l 

was estimated in the previous Lemma, and 

the following upper bound: 

t N \ 



233(^,{a)) n*^ 

Vj=i ) 

was established for it. The second difference can be expanded as: 

( N N \ N (( 



e 2 



N 



\ 7=1 7=1 / k=\ W ]<k j 



*(Pfk-fk)* 



V j>k j 



SO we can write: 



( N 



N \ 



k=l j<k j>k 



yj=i 7=1 J 

and, using Theorem |7.9[ we can establish that this is less than or equal to 



I N 



k=\ j^k J\j=l 



e V2 



N 



(264) 



(265) 



(266) 



(267) 



yV7=i / 



If Ck(^) = supyy A/''"'''^'e ^ then the absolute value of the difference of Eq. 265 is less than or equal than: 

C.,.,(^) ^^^ ^ 

l^\aj\+l 

\k=l j^k J V 7=1 J 



( N 



(268) 



and the statement follows. 



8 The noncommutative Kubo formula on the discrete torus 

The noncommutative Kubo formula with a dissipation F on the torus is a straightforward translation of the exact Kubo 



formula from Eq. 87 This translation can be achieved using the maps P and / defined in Eqs. 197 and 198 and the 



approximate differential calculus introduced in the previous section. It takes the form: 



d-ij{N) = Tt (m) * (P o r o / + Xa)-'5,<I)pd(^)) , 



(269) 



where h = (P o ^^/^^)h, h being the Hamiltonian in the thermodynamic limit. We aie now in the position to make two 



important statement about the above approximate Kubo. The first statement is that Eq. 269 can be efficiently evaluated 



on a computer, to a point where, for example, it enabled us to gain unprecedented insight into the transport properties 
of the Quantum Hall systems and topological insulators. The second statement is that, at least for the relaxation time 
approximation F = l/r,-, the errors decay exponentially fast as the size of the torus is increased. As such, the bulk 
transport coefficients can be computed with good accuracy on relatively small simulation boxes. 
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8.1 Error bounds 



We will assume the relaxation time approximation, where F is equal to l/r,- times the identity map. The parameter 
T,. < oo is the so called relaxation time, an empirical parameter that is usually extracted from the experimental data. 

Theorem 8.1. Let k be defined as: k = min{l/2r, , TikT} - e, with e an arbitrarily small but strictly positive constant. 
Then, for the nearest-neighbor hopping Hamiltonian ofEq. [72] there exists a finite and fully identifiable constant 93(^) 
such that: 

|cr„-^,/yV)|<23©6'-5f^, (270) 

for all < ^ < sinh-\K/2d). 

Proof. The statement will follow after we put together the comparative estimates derived in the previous sections. 
We will follow our general strategy and compare the exact noncommutative Kubo formula with an approximate Kubo 
formula from the periodic algebra, which will be subsequently compared with the approximate Kubo formula from the 
torus algebra. We start from the exact conductivity tensor: 



and use Proposition 4.2 



aj, = T({djh) * (l/r, + Xft)-'5,OFD(/t)) , 
in reverse to write: 



cr 



jk 



Xoo 
dte-'l-'T({djh) * M„(f,0)[5,(DFD(A)]) . 



With the explicit representation of the time evolution given in Eq. 46 the above expression becomes: 

dte-'l''T({djh) * e"'"' * [d^^^vTiQi)] * e'"') 



(271) 



(272) 



(273) 



Now one can see that we can connect with the Theorem 6.3 which allows us to write the following comparative 
estimate on the integrand: 



It (idjh) * * [dkMm * e'"') - Tper (idjh,,) * e-"'" * [5,OFD(ftp)] * e"'")\ 



( 4 



<23i(^, [a]) 



N-'e-'f^''. 



(274) 



where hp = ^\^h. The O's can be easily identified as: <l>i(z) = z, ^i^iz) = e"^'^' and O3 = Ofd(z)- 

Furthermore, with h defined as h = (P o S\i^^)h, one can see that hp = Jh. Hence, we can use Theorem 
write: 

kper ((djhp) * e-'"^ * [d,(^PD{hp)] * e""") - Tt ({djh) * * [5,(Dfd(/i)] * e"H 



7.12 



to 



< 934(^,{a)) 



(4 A 



(275) 



where O/s are same functions as before. 

With K chosen as indicated in the text of the Theorem, the complex set defined by those z's for which 
dist(z, cr( Jh)) < k (which is a subset of the complex strip |lm(z)| < k) resides inside the analytic domain of the functions 
4>,(z)'s. In fact the only function that needs to be checked is O3, the Fermi-Dirac function, which has poles at 
Ep + i(2n + \)nkT, « e Z. The functions 02,4(2) depend on the parameter f which is to be integrated from to 00, 
with a weight factor e"'^""' . Since, 02,4(0 ^ const, e"', the integration over f remains finite when k is chosen less than 
1 /2Tr. We then conclude: 



dt e-'I'^T ({djh) * e''"' * [40fd(/!)] * - J dt e-'>^'Tf ((djh) * e"'"' * [5iOFD(^)] * e'''') 



const. 



1 /Tr - 2k 

and the statement follows. 



(276) 
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Hall plateau 




Pl2 



Magnetic Field 
Figure 5: Typical output of a Hall measurement. 

8.2 Numerical Implementation 

The only issue remaining to be clarified is how to invert the map 1 /r, + JLj^ will restrict the discussion to the 
relaxation time approximation). Formally, this does not present any fundamental difficulty because we can simply 
assume to work in the Hilbert space defined by the scalar product: 



(L~g) = rr(r*~g), nil) 

which was already discussed in our introductory sections. It is straightforward to demonstrate that X;, becomes a 
bounded normal operator (more precisely, /X^ is self-adjoint) so the inverse (1/t + X/,)"' can be computed using the 
ordinary spectral decomposition of X/v This is precisely the route we chose to follow here. In the following, we explain 
how to compute the spectral decomposition of X/v 

We first introduce an equivalence relation on O/,^ by saying that w ~ w' if there is a g e such that co = r'^co' . 
We will denote the equivalence classes by [oj\. The factorization of Q^v in equivalence classes will be denoted by 
Q.nI ~, and the factorization of the probability measure dPf(to) by dP-rlaj] (= dPf{oj)l ~). 

Using Proposition |7. 3 1 the scalar product can be written as: 

(/>§) = ^ r dPr{u)lxUnJ*){nM, (278) 

and due to the covariant property of the representation n, we observe that: 

Tr|(7r„f)(7r^|)) (279) 

remains unchanged when a> is replaced by XgCj, with q an arbitrary point of the torus. In other words, the trace in Eq. 90 
depends only on the equivalence class of oj so we can conclude: 

(/,§)=^r dPT[aj]Ti{{nJ*)i7T^~g)], (280) 

where a> appearing inside the trace can be any element from the equivalence class [oj]. Furthermore, since the tt^/ 
representation of an element involves only the values of / on [w] x Tj, we can draw a more fundamental conclusion, 
namely, that the original Hilbert space decomposes as a direct integral: 



'H= ^/PtM-^m (281) 
where the elements of the Hilbert space are functions / : [oj\ x — > C and the scalar product is: 



(/,§)[„] =Tr{(7r„r)(ir^|)|. (282) 

Here, again, the o) appearing inside the trace can be any element from the equivalence class [cj]. The product and the 
star operations of the torus algebra remain well defined when restricted to "//[i^]. 
Now, the operator X/, leaves invariant the Hilbert spaces so 



£f,= dPrlu] £r,[ojl (283) 
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where £.},[lo] is simply the restriction of the map X/, to the space of functions defined over [co] x Td- The spectral 
decomposition of can be achieved in the following way. We pick a representative co from the class [oj] and 

instruct a computer to diagonalize the operator n^^h. Since the Hilbert space ^^(T^) is finite, this can be achieved with 
the diagonalization routines from any standard linear algebra library. Let {e^, 0") be the resulting eigensystem: 

(7fJ)C = e«C> a = l'---'|Trfl- (284) 

Now note that all oj's from [co] can be uniquely written as XgO), with q a point of the torus. For p and q arbitrary points 

of the torus, we define: 

fab(j,u,p) = raiq)W^y"^"''^'"'^^ a,b=l,..., ITrfl. (285) 



Eq. 285 defines genuine functions from [to] x Tj to C, which will be viewed as elements of ^iioj]. A straightforward 
calculation will show that: 

1 . The elements are orthonormal: 

(fah,fcd)[aj]=SacSbd- (286) 

2. The sequence {fah}a,b=i,...,\T,,\ is complete, for if: 

if,fab)M=0 (287) 

for all fl, ft = 1, . . . , |Trf|, then / is necessary the null element. 

3. The elements are eigenvectors of X/,[w]: 

£d<^](U) = i{ea - eb)fab- (288) 

4. The operator representation of {fab\a,b=\,...,\Tj\ is: 

TiJah = ICX-A^I (289) 
Thus, we have accomplished the spectral decomposition of and we can conclude: 

(1/T, + X,M)-'(/) = V (290) 
^ I t,. + i{ea - eb) 

a,b=\ 

Given that n^^fab = 10a we can now write the explicit expression for the conductivity tensor: 

r r n 1 V <CI^-(^;^)IC)<'A^I^«(4%d(^)IC) 

0-jk = - dPjloj] — > — . (291) 

J Trf ^ I T,. + i{ea - eb) 

This exact formula was implemented on a computer and used in our numerical simulations. 

We should comment here that computing the inverse (l/r^ + -Cfjlco])'^ is, after all, a classic linear algebra problem 
so it can be achieved in many different ways. All that is required is a clever implementation of the action of on 
the elements of We are currently exploring various other paths to invert 1/t,. + £.i,[a)], but for the present work 
we chose to use the spectral decomposition because this way we can understand and minimize the numerical errors. In 
fact, the computations of o-jkiN) shown in the next section are done with the machine precision. 



9 Application: Tlie Integer Quantum Hall Effect 

There is a tremendous amount of literature dedicated to the Integer Quantum Hall EflFect (IQHE), so here we will 
be brisk and stay to the point. The reader can find a comprehensive discussion of the effect in the classic textbook 
put together by Prange and Girvin [46 1. Most of the facts stated below can be found in this textbook. There have 
been several important developments after the publication of this book, namely the observation of IQHE in graphene 
ir70ll4ni42l|43ll , in 2-dimensional Quantum spin-Hall systems [29 1 and at the surface of a 3-dimensional topological 
insulator [fT3ll22ll65l . It is particularly important to mention that IQHE in graphene was observed at room temperature 
[|43l . hence making our finite temperature analysis even more relevant. An explicit lattice analysis of the IQHE in 
clean graphene can be found in Refs. ifTTl iTl and the effect of disorder (and interaction) at zero temperature has been 
analyzed in [58 1. Our goal here is to demonstrate what we can achieve with the new numerical procedure. The analysis 
of the IQHE itself, based on the new numerical results, is beyond the scope of this work, simply because that by itself 
could be an entirely separate research project. 
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Figure 6: The density of states of our model Hamiltonian as function of energy and magnetic flux through the unit cell, 
for three degrees of disorder: W=0, 2 and 3. At W=0, one can recognize the Hofstadter butterfly. 




-4 o 4-4 o 

Energy Energy 



Figure 7: DOS (left) and color map of the Hall conductivity (right) for W = 0. The regions of quantized Hall 
conductivity, which appear as well defined patches of same color and, are indicated at the right. 



In short, IQHE is observed in the transport measurements on a 2-dimensional electron gas when subjected to a 
perpendicular magnetic field. Experimentally, the 2-dimensional electron gas can be obtained by trapping free-moving 
electrons at the interface between two semiconducting materials. A 2-dimensional electron gas is also forming at the 
metallic surface of a 3-dimensional topological insulator and it occurs naturally in the atomically thin graphene. The 
density of the 2-dimensional electron gas can be adjusted in a finely controlled manner by applying gate voltages. If 
the samples are of high quality and the temperature is low enough (except for graphene where IQHE can be observed at 
room temperatures in very high magnetic fields), the transport curves of the Hall resistivity pi2 (to be defined shortly) as 
function of the B-field (while holding the electron density constant), or as function of electron density (while holding B 
constant) display a staircase pattern as schematically illustrated in Fig. |5] with the horizontal plateaus being quantized 
to several digits of precision. In the same time, the diagonal resistivity pn is observed to vanish inside the Hall plateaus, 
and to have sharp peaks at the transitions between plateaus. The experimental data for IQHE is most often presented 
in terms of the resistivity tensor p,y = (cr"^),j, rather than the conductivity tensor o-jj, because the resistivity can be 
measured with much greater precision. Explicitly, for our isotropic model; 

Pi I = 0"! 1/(0-2 +0-2 ) 

(292) 

P12 = 0-12/(0-11 + 0-12). 

There is a special feature in 2-dimensions, in that the resistivity has same physical units as the resistance and 
conductivity as conductance. 

In the absence of a periodic potential, the energy spectrum of the 2-dimensional electron gas (in a mean-filed, 
non-interacting approximation) consists of equally spaced Landau levels separated by spectral gaps. One can easily 
compute the Hall conductance an at T = 0, assuming the lowest n Landau levels as completely filled, in which case the 

2 

result will be an = rij = nGo. While this can provide a clue about the origin of the quantized values, it cannot explain 
the existence of the plateaus. That is because the probability of having some Landau levels completely filled while the 
rest completely empty is practically zero. There is no way one can adjust the electron density with such precision! Also, 
within such a simplistic picture, it is impossible to explain the variation of the electron density when traversing a Hall 
plateau. Presently, the accepted view is that, in between the Landau levels, there are localized impurity states, in which 
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Figure 8: DOS (left) and color map of the Hall conductivity (right) for W = 2. The regions of quantized Hall 
conductivity, which appear as well defined patches of same color and, are indicated at the right. 



case the Fermi level transitions smoothly (instead of jumping) from one Landau level to another, when increasing the 
B-field or the electron density. The locahzed spectrum is thought to be generated by a disordered potential, inherently 
present even in high quality samples. Thus, the presence of disorder is crucial for the explanation of the IQHE. A 
complete, rigorous theory of the IQHE at zero temperature was derived in a series of papers by Bellissard and his 
collaborators. In particular, Ref. [6| developed the noncommutative theory of transport that we followed here, together 
with a set of optimal conditions for the quantization and invariance of the Hall resistance in the presence of strong 
disorder. 

The hypothesis of localized states was experimentally tested in several ways, and the conclusions were always 
positive. For example, it was observed that the width of the plateaus increases when the mobility of the electrons at 
zero field decreases. However, too much disorder can lead to the disappearance of the plateaus. Other interesting and 
relevant experimental observations is that the shape of the curves at the beginning or at the end of a Hall plateau is not 
always smooth as illustrated in Fig. |5] Dips and sometime even oscillations can show up in the experimental curves. 
There is also a very interesting dependence on temperature of the transport coefficients. One standard observation is the 
gradual appearance and sharpening of the Hall plateaus as the temperature is lowered. The direct resistance displays 
an activation behavior which can be understood from the simple shape of the Fermi-Dirac distribution, but when the 
temperature becomes extremely low a different behavior sets in, determined by the type of disorder and dissipation 
present in the system. 

There are many published simulations of the IQHE in the presence of disorder ifTSl [161 [36l l63l l56l ISTl l68l l57l 
l62l ItTI l38l l32l l58l l35l [TTl . However, many of these simulations are restricted to just the diagonal component of the 
conductivity tensor. The computation of o-\2 or pn is more difficult and, as such, simulations of (T\2 or are scarce. 
Almost all of the existing simulations present o-\2 as function of Fermi energy at fixed magnetic field. Such plots, 
of course, show wide and well quantized Hall plateaus when the Fermi energy crosses the gaps between the Landau 
levels. In experiments, however, what is controlled is the electron density and not the Fermi level. Since the variation 
of the electron density is very small when the Fermi level crosses the gap between the Landau levels, the wide and well 
quantized Hall plateaus would appear much narrower if an was plotted as function of electron density. In fact, it is not 
clear at all if any of the Hall plateaus observed in these numerical studies would remain visible if the data was plotted 
as function of electron density. Furthermore, we were not able to find any record of a simulation which analyzes the 
conductivity or resistivity tensor as function of magnetic field at fixed electron density. 

We now present what we were able to achieve (so far) with the formalism presented in the previous sections for 



the model Hamiltonian of Eq. 12 As we already mentioned, at the algebra level, this Hamiltonian is generated by h 



which takes the values h(oj, m) = 1 if |m| = 1, h{oj, 0) = VKwo, and zero in rest. The calculations will be carried out on 
a finite lattice which is wrapped into a torus. The size of the finite lattice will be indicated as Nr x Nr, where Nr is 
the total number of nodes in one direction. In the old notation, Nr would be given by IN + 1 but here we relax the 
requirement that Nr be an odd number. In two dimensions {d = 2), the matrix F is simply Fij = ei^cf), where cp is the 
magnetic flux through the repeating cell, expressed in units of flux quantum 0o = h/e. The Hamiltonians at fluxes 
and + 1 are unitarily equivalent so we can restrict the values of cp to the interval [0, 1]. Since the calculations are 
carried out on the torus, the magnetic flux will always be quantized in our calculations according to the rule given in 



Eq. 182 As such, the larger the torus the more magnetic field values we can sample. 

We start with some qualitative analysis. The spectrum of h, obtained on a 120 x 120 lattice, is shown in Fig.|6]for 
three values of W: 0, 2 and 3. The spectrum is represented via the smoothed density of states (DOS), defined as: 



D(e) = ^T(lm(h-e-i5r'), 



(293) 
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Table 1: The numerical values of crn at kT = l/Tr = 0.1, W = ^ = and various Fermi energies, obtained with the 
new algorithm for increasing lattice sizes. The last column displays the value of crn computed with machine precision 
via Eg. [2941 



Ef 


80 X 80 


100 X 100 


120 X 120 


140 X 140 


Exact 




4.0339628247 


4.0339630615 


4.0339630708 


4.0339630712 


4.0339630712 


-0.4 


3.9394154619 


3.9394154735 


3.9394154621 


3.9394154624 


3.9394154624 


-0.8 


3.7040304262 


3.7040301193 


3.7040301310 


3.7040301307 


3.7040301307 


-1.3 


3.3684805414 


3.3684801617 


3.3684801517 


3.3684801516 


3.3684801516 


-1.7 


2.9522720814 


2.9522713926 


2.9522714007 


2.9522714009 


2.9522714009 


-2.2 


2.4678006935 


2.4678005269 


2.4678005093 


2.4678005104 


2.4678005104 


-2.6 


1.9239335953 


1.9239338070 


1.9239338090 


1.9239338089 


1.9239338089 


-3.1 


1.3274333126 


1.3274333067 


1.3274333084 


1.3274333085 


1.3274333085 


-3.5 


0.6854442914 


0.6854442923 


0.6854442923 


0.6854442923 


0.6854442923 


-4.0 


0.1086465150 


0.1086465150 


0.1086465150 


0.1086465150 


0.1086465150 



Table 2: The numerical values of crn at kT = 1/t,. = 0.025, W = ^ = and various Fermi energies, obtained 
with the new algorithm for increasing lattice sizes. The last column displays the value of crn computed with 



machine precision via Eq. 294 The numbers in parentheses represent the relative errors, computed as Ictj i (approx) 
o"! 1 (exact)l /cTi 1 (exact). 



Ef 


100 X 100 


120 X 120 


140 X 140 


160 X 160 


Exact 


0.0 

.4/, 
-2 X 1/9 

-3 X ♦/■) 

-4 X 1/9 
-5 X 4/9 
-6 X 1/9 


16.204612406 (3.7e-5) 
15.800901211 (8.2e-5) 
14.845640534 (8.3e-5) 
13.501282967 (1.7e-4) 
1 1.830160006 (l.le-5) 
9.8929640489 (3.3e-4) 
7.7119398206 (2.4e-5) 


16.204972260 (1.5e-5) 
15.799191631 (2.5e-5) 
14.847746111 (5.8e-5) 
13.498925093 (1.3e-6) 
11.830983739 (5.8e-5) 
9.8886679730 (l.Oe-4) 
7.7130531856 (l.le-4) 


16.205119575 (6.6e-6) 
15.799439976 (9.7e-6) 
14.846643808 (1.5e-5) 
13.498688837 (1.6e-5) 
11.830916078 (5.2e-5) 
9.8899094634 (2.2e-5) 
7.7127298903 (7.7e-5) 


16.205181208 (2.8e-6) 
15.799694877 (6.3e-6) 
14.846886090(6.16-7) 
13.498843113 (4.6e-6) 
11.830155334 (1.16-5) 
9.8896881580 (4.7e-7) 
7.7123144240(2.36-5) 


16.205227112 
15.799593904 
14.846876954 
13.498906270 
11.830294392 
9.8896834893 
7.7121326232 


-7 X */9 
-8 X 1/9 
-9 X 4/9 


5.3242604491 (2.4e-5) 
2.7469871332 (2.4e-4) 
0.1099066152 (2.8e-9) 


5.3236541115 (1.3e-4) 
2.7475442292 (4.0e-5) 
0.1099066156 (Oe-10) 


5.3241123103 (5.2e-5) 
2.7476748910 (7.4e-6) 
0.1099066156 (Oe-10) 


5.3243169737 (1.3e-5) 
2.7476680440 (4.96-6) 
0.1099066156 (Oe-10) 


5.3243912832 
2.7476545340 
0.1099066156 



Table 3: The numerical values of crn at kT = I/t^ = 0.01, W = (p = and various Fermi energies, obtained with 
the new algorithm for increasing lattice sizes. The last column displays the value of crn computed with machine 
precision via Eq. |294[ The numbers in parentheses represent the relative errors, computed as |crn(approx) - 



crn (exact)|/crn (exact). 



Ef 


100 X 100 


120 X 120 


140 X 140 


160 X 160 


Exact 


0.0 

-1 X 4/9 
-2 X 1/9 
-3 X 4/9 
-4 X 1/9 

-5 X 4/9 
-6 X 4/9 
-7 X 4/9 

^-8x1/9 

-9 X 4/9 


40.520686410(1.26-4) 
39.517833421 (3.3e-4) 
37.110642808 (2.96-4) 
33.777005514 (7.7e-4) 
29.585203369(2.16-4) 
24.787127553 (2.4e-3) 
19.265355553 (8.96-4) 
13.320411822 (5.5e-4) 
6.8441390250 (3.96-3) 
0.1101512994 (2.3e-5) 


40.522558677 (7.56-5) 
39.497150004(1.96-4) 
37.138820375 (4.66-4) 
33.748897618 (5.4e-5) 
29.596597708 (6.06-4) 
24.713893823 (5.2e-4) 
19.304253858 (1.16-3) 
13.288432535 (1.8e-3) 
6.8590106166(1.76-3) 
0.1101540276 (1.7e-6) 


40.523641888 (4.86-5) 
39.499629103 (1.2e-4) 
37.114579007 (1.86-4) 
33.742729517 (2.3e-4) 
29.596944121 (6.16-4) 
24.729175598 (9.6e-5) 
19.305176839(1.16-3) 
13.302525433 (7.8e-4) 
6.8736981172 (3.96-4) 
0.1101538290 (7.9e-8) 


40.524308038 (3.26-5) 
39.508996214(1.06-4) 
37.120997917 (1.26-5) 
33.747792038 (8.7e-5) 
29.573808787 (1.66-4) 
24.724104241 (l.Oe-4) 
19.288293494 (2.96-4) 
13.305757131 (5.4e-4) 
6.8765917049 (8.16-4) 
0.1101538378 (Oe-10) 


40.525626855 
39.504691905 
37.121462109 
33.750748672 
29.578697042 
24.726791044 
19.282623910 
13.313040334 
6.8710107299 
0.1101538378 



where 6 was given the small value 0.01. The DOS is plotted as an intensity map in the {e,(f>) plane. The plot atW = 
reproduces the famous Hofstadter butterfly [23 |. This can serve as a reminder that the simple Landau level picture 
becomes more complicated when a periodic potential is present. The Hofstadter spectrum has a fractal structure with 
the patterns seen in the first panel of Fig. |6]repeating over and over when zooming to finer and finer scales. At weak 
magnetic fields, and especially towards the edges of the energy spectrum, one can identify well defined Landau bands 
(a more careful scrutiny will reveal that these bands are made of many sub-bands). If the Fermi level is somehow 
fixed in between these bands, then one expects quantized values of cri2 in units of Go- The DOS atW = 2 and W = 3 
look like a blurred version of the DOS atW = 0. The blur originates from the impurity states generated inside the 
empty spaces between the Landau bands by the disordered potential. Maps of the Hall conductivity are reported in 
Figs.|7]|8]and|9]for W = 0,2, and 3 cases, respectively. These calculations were performed on a 100 x 100 lattice with 
kT = 1/t,- = 0.01. A single disorder configuration was used, since the fluctuations due to the disorder will not have 
any visible eflFects for the intensity maps of Figs.|6]|7]|8]and|9j The Fermi energy was varied over the entire energy 
spectrum, which was sampled at 60 equally spaced points. In all three cases, the maps display regions where cri2 takes 
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Figure 9: DOS (left) and color map of the Hall conductivity (right) for W = 3. The regions of quantized Hall 
conductivity, which appear as well defined patches of same color and, are indicated at the right. 
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Figure 10: Diagonal and Hall resistivities computed with Eq. 295 for 10 random disorder configurations, at disorder 
strengths W = I, 2 and 3. The simulation was performed on a 100 x 100 lattice with kT = I /t^ = O.l. The electron 
density was fixed at n^, = 0.25. The standard deviation was observed to be less than 2%. 



the expected quantized values. In the middle of these regions, the quantization occurs with five digits of precision or 
better. For the cases with disorder, it is important to notice that this quantization occurs beyond the regions empty 
of spectrum, that is, in the regions where Ef is emersed in the localized spectrum. The Hall conductivity maps from 
Figs. [7] |8] and |9] give a useful panoramic view of the Hall conductivity but, even though we see sharp quantized values, 
the results are actually far from demonstrating the IQHE in our model. That is because in the real experiments, there 
is no control over the Fermi level. What is being controlled is the density of electrons and the Fermi level adjusts itself 
correspondingly. What it is known for sure is that the Fermi level will never stay in the regions empty of spectrum. As 
such, to actually demonstrate the occurrence of the IQHE in our model, we need to map the transport coefficients as 
function of density while holding the magnetic flux fixed, or as function of magnetic flux while holding the electron 
density fixed, allowing the Fermi level to adjust itself to the appropriate value. Revealing the Hall plateaus in these 
new conditions is a much more difficult task. At last, it is interesting to note the similarity of the Hall conductivity map 
at W = in Fig. [7] and the Chern number map computed in Ref [44 J. There are expected difl'erences between the two 
results because the map in Fig. |7] was computed at a finite temperature while the Chem number map gives the Hall 
conductivity at zero temperature. 

At this point we start our quantitative analysis. First, we present a simple but direct test of the numerical algorithm. 
In the absence of disorder and magnetic field, the Kubo formula takes the extremely simple form: 

(Tn = :r-p^ dh dk2 (294) 

where Ek is the band energy of the model: E^ = 2 cos ki +2 cos ^2- This formula can be computed with the machine 
precision and we will call such a computation "exact." In Tables [T] [2] and |3] we report a comparison between the 
exact value of crn and the values obtained with the new algorithm on lattices of varying sizes. In Table [T] we fixed 
the temperature and dissipation at relatively large values: kT = 1/t,. = 0.1. Four lattice sizes were considered: 80 x 80, 
100 X 100, 120 X 120 and 140 x 140. The Fermi level was varied from the bottom to the middle of the energy spectrum. 
By examining the columns of Table [T] we see that, already for the smallest lattice size, the output from the new 
algorithm displays 7 digits of precision for all E/s. As predicted, the degree of precisions increases rapidly with 
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Figure 11: First row (Second row): The diagonal and the Hall resistivities as function of Fermi energy (density) at 
fixed magnetic flux (p = O.lcpo, for kT = I/t^ = 0.025 and disorder strengths W = 1, 2, 3. Each panel compares the 
data obtained on the 100 x 100 lattice (circles) and on the 120 x 120 lattice (squares). 



the lattice size and for the largest size we see 11 digits of precision for all Ef's. This indicates that all the lattice 
sizes appearing in Table [1] are in the asymptotic regime where the exponentially-fast convergence occurs. In Table [2] 
the temperature and dissipation were fixed at lower values: kT = l/r^ = 0.025. Four lattice sizes were considered: 
100 X 100, 120 X 120, 140 x 140 and 160 x 160. In this table too, we see an improvement of the accuracy as the size 
of the system is increased, though not as dramatic as in Table [T] Looking at the relative errors, we see that, on average 
(with respect to Ep), the accuracy was improved by more than one order of magnitude when the lattice size was 
increased from 100 x 100 to 160 x 160. For the largest lattice size, there are more than 5 digits of precision. Based on 
the monotonic reduction of the relative errors from one column to the other (with one or two exceptions), we believe 
that the asymptotic regime where the exponentially fast convergence occurs has been already reached in Table |2] In 
Table|3] the temperature and dissipation were fixed at even lower values: kT = I/t^ = 0.01. Here we still see a weak 
improvement of the accuracy with the size of the system. The precision here is about 4 digits for the largest two lattices. 
Based on the information contained in these tables, we expect the algorithm to be well converged on the 120 x 120 
or 140 X 140 lattices when kT and 1/t,. are equal to or larger than 0.025, and we expect reasonable outputs even for 
kT = 1/t,- = 0.01. We will always verify these statements explicitly, by repeating the calculations for increasing lattice 
sizes. 



Next we discuss the issue of averaging over disorder. The integrand in Eq. 291 



1 



|T,,| 

z 

a,b=l 



(0^|^^(g/i)IO<.A^|7r^(a(DFD(^)|<^P 
1/t + /(6„ - efo) 



(295) 



is self-averaging, meaning it converges to its disorder-average as the size of the torus is taken larger and larger. This 
is an important observation because it practically solves the issue of disorder average. For the sizes considered in this 
study, we can demonstrate that indeed, the integrand is practically independent of the particular disorder configuration. 



so we can ehminate the disorder-average step entirely. Fig. 10 presents a computation of pn and of pi2 as functions of 
the magnetic flux for 10 randomly generated disorder configurations, and for three disorder strengths, W = 1,2 and 3. 
The calculation was performed on a 100 x 100 lattice at fixed electron density He = 0.25, and at kT = l/r,- = 0.1. The 



standard deviation of the fluctuations due to disorder is less than 2% for the data reported in Fig. 10 This calculation 
confirms that, starting from this lattice size, the fluctuations of the resistivities from one disorder configuration to 
another are very small, at least for this temperature and dissipation strength. From now on, the transport coefficients 



will be computed with just one disorder configuration, using the formula from Eq. 295 and we will visually estimate 
the fluctuations due to the disorder by examining the data for diff'erent lattice sizes. 
We now present the actual simulation of the IQHE. The first row of Fig. 



11 reports pn and pi2 as function of 



Fermi energy, with the magnetic flux hold fixed at O.lcpo and for kT = I /r,. = 0.025. The calculations were performed 
on the 100 X 100 lattice and then repeated on the 120 x 120 lattice, for three disorder strengths: W = 1,2 and 3. The 
most important observation here is that the data corresponding to the two lattice sizes overlap almost perfectly for all 
three disorder strengths, indicating a very good convergence of the results. Also, there are no visible fluctuations due 
to disorder. For = 1 and = 2 we can see four well defined and quantized Hall plateaus. The diagonal resistivity 
takes very small values in the region Hall plateau regions. For W = 3, only the first two Hall plateaus display good 



quantization. The same data on the resistivity tensor is plotted in the second row of Fig. 1 1 as function of density 
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Figure 12: First row (Second row): The diagonal and the Hall resistivities as fvmction of Fermi energy (density) at fixed 
magnetic flux ^ = O.150o, for kT = I/t^ = 0.025 and disorder strengths W = 1,2, 3. For clarity, we show only the data 
obtained on the 120 x 120 lattice. 



W = l W=2 W = 3 
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Figure 13: The diagonal and the Hall resistivities as function of magnetic flux, at kT = l/r,. = 0.1, = 0.25 and 
disorder strengths W = 1,2, 3. The first row shows a comparison between the data obtained on the 100 x 100 lattice 
(circles) and on the 120 x 120 lattice (squares). The two sets of data shown in the first row were joined together and 
plotted in the second row. 



He =T (l/(l + e^^'^''"'!'^^)^. Note that rig varies from (no electrons) to 1 (one electron per lattice site). As we already 
pointed out, the quantized Hall plateaus are drastically reduced in this representation. The disorder seems to slightly 
widen the Hall plateaus, especially when one compares V7 = 1 and W = 2 plots, but we do not see the anticipated 
dramatic eff'ects of the disorder. Calculations on even larger lattices (see the following discussion) makes us confident 
that the calculations shown in Fig. 1 1 are converged. Additional calculations will also show that there are no visible 
changes as the temperature and dissipation are lowered, so the narrow plateaus seen in our simulations must be a 
characteristic of the model at the particular value of the magnetic flux considered here. Looking back at Figs. [8] and 
|9] one can see an increase of the spreading of the Landau bands with the disorder when the magnetic flux approaches 
0.5^0' so we have repeated the calculations by fixing the magnetic flux at the larger value (p = O.150o. in a hope to 
improve outcome. The results are shown in Fig. 12 and, indeed, here we can see the first and second quantized Hall 
plateaus forming when the disorder is increased, even when the data is plotted as function of electron density. Still, 
we think there is a lot of work to be done until we will see a quantitative similarity with the experimental data. To our 
knowledge, this has been never achieved. 

Next, we fix the electron density and vary the magnetic flux. Fig. 13 reports pn and pu as functions of the 
magnetic flux, with the electron density hold fixed at 0.25 and for kT = 1/t,. = 0.1. The calculations were performed 
on the 100 X 100 lattice and then it was repeated on the 120 x 120 lattice, for three disorder strengths: W =1,2 and 3. 
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Figure 14: The diagonal and the Hall resistivities as function of magnetic flux, at kT = 1/t,. = 0.025, rig = 0.25 and 
disorder strengths W = I, 2, 3. The first row shows a comparison between the data obtained on the 100 x 100 lattice 
(circles) and on the 120 x 120 lattice (squares). The second row shows a comparison between the data obtained on the 
120 X 120 lattice (squares) and on the 140 x 140 lattice (triangles). The two sets of data shown in the second row were 
joined together and plotted again in the third row. 



The results are compared in the first row of Fig. 13 We see again a very good overlap between the data corresponding 
to the two lattice sizes, for all three disorder strengths, and there are no visible fluctuations due to disorder. For each W 
values, we can then gather together the data obtained on the two lattices together and practically double the magnetic 



flux sampling values. The resulting data is plotted in the second row of Fig. 13 and we believe these graphs are quite 
accurate representation of the resistivity tensor of the model system, at this temperature and degree of dissipation. The 
features seen in these graphs, especially those for W = 2 and 3, resemble qualitatively the experimental IQHE curves 
at higher temperatures, with the Hall plateaus just about to form and with the diagonal resistivity dipping in the plateau 
regions. 



In Fig. 14 the temperature and dissipation were reduced to kT = I/t^ = 0.025. All the other parameters were 
held at the same values as above. The first row presents a comparison between the data obtained on the 100 x 100 
and 120 x 120 lattices, and the second row presents a comparison between the data obtained on the 120 x 120 and 
140 X 140 lattices. The overlap between the two sets of data is already very good in the first row, and it is almost perfect 
in the second row. The visible discrepancies are due to different sampling of the magnetic flux. We join together the 
data from the 120 x 120 and 140 x 140 lattices in order to increase the magnetic flux sampling values, and the plot of 



the resulting data is shown in the third row of Fig 14 Examining the Hall resistivity, we see the first Hall plateau to be 
very well defined and relatively flat, at all three disorder strengths. The first plateau is seen to widen out as the disorder 
strength is increased. The calculations also resolved higher order Hall plateaus. For W = 1 and 2, we distinguish as 
many as five Hall plateaus. For all these five Hall plateaus (excepting the 4th one), the diagonal resistivity dips to very 
low values, as expected. The fifth and fourth Hall plateaus are washed out when disorder is increased toW = 3. 



Finally, in Fig. 15 the temperature and the dissipation were reduced even more to kT = 1/t,. = 0.01. The layout 



of this figure is the same as for Fig. 14 and there are practically no major changes from the previous case, except that 
the first Hall plateau became sharper. In fact, a direct comparison shows a good quantitative similarity between the 
data at kT = 1/t,. = 0.025 and kT = I/t^ = 0.01, when the graphs were visually inspected. 
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Figure 15: The diagonal and the Hall resistivities as function of magnetic flux, at kT = I/t^ = 0.001 and disorder 
strengths W = 1, 2, 3. The first row shows a comparison between the data obtained on the 100 x 100 lattice (circles) and 
on the 120 x 120 lattice (squares). The second row shows a comparison between the data obtained on the 120 x 120 
lattice (squares) and on the 140 x 140 lattice (triangles). The two sets of data shown in the second row were joined 
together and plotted again in the third row. 

10 Conclusions 

The framework of the C'-algebras and noncommutative calculus enabled us to develop a "canonical" finite-volume 
approximation to the exact noncommutative Kubo formula with dissipation. This approximate formula is amenable 
on a computer. Within the same framework and for the relaxation time approximation for the dissipation effects, we 
were able to establish rigorous error bounds for this approximate formula, and to demonstrate that the errors vanish 
exponentially fast in the thermodynamic limit. Given the general upper bounds derived here for generic correlation 
functions, we believe that this last conclusion also applies to situations beyond the relaxation time approximation. The 
C* -algebra formalism was absolutely instrumental for both aspects of our analysis, namely, for the derivation of the 
approximate Kubo formula and for the analysis of the errors. We believe that the formalism of C'-algebras will also 
lead to better numerical implementations of the finite-volume approximate Kubo formula. 

Based on the abstract theoretical ideas, we devised a numerical algorithm for the computation of the conductivity 
tensor for disordered systems under magnetic fields. The main computational advantage of the algorithm is that it does 
not require twisted boundary conditions and integration over a Brillouin torus in order to achieve exponentially fast 
convergence for the thermodynamic limit. The time saved by not having to repeat the calculations for many diff'erent 
twisted boundary conditions was fully re-invested to increase the size of system, leading to a substantial boost of the 
simulation box over the existing simulations. 

We presented an application to the Integer Quantum Hall Eff'ect. Working with a lattice model of a disordered 
2-dimensional electron gas in perpendicular magnetic fields, we were able to compute the resistivity tensor on lattices 
containing as many as 140 x 140 sites, placing the calculations quite close to their thermodynamic limit for the 
temperatures and the degree of dissipation considered in our study. We computed the resistivity tensor while varying 
the Fermi energy (or equivalently the electron density) while holding the magnetic field constant, and as function of 
the magnetic flux while holding the electron density fixed. We were able to resolve as many as five Hall plateaus, 
with the first plateau appearing sharp and well quantized as the temperature and dissipation is reduced. We could also 
witness how the first and second plateaus become wider and wider as the disorder strength is was being increased. To 
our knowledge, this is the first time such simulation results have been reported. 
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